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SUMMARY 


This  report  describes  the  results  of  two  distinct  studies  of  the 
aerodynamic  heating  of  nose  tips.  In  PART  1  of  these  studies  we  examined 
experimentally  the  development  of  laminar,  transitional  and  turbulent 
boundary  layers  over  ablated  nose  shapes  at  Mach  numbers  from  8  to  13  for 
Reynolds  numbers  up  to  80  x  10^  based  on  model  diameter.  Two  groups  of  models 
were  employed  in  this  study,  blunt  elliptic,  blunt  biconic  and  triconic  models 
were  used  in  the  first  phase  of  this  investigation  where  the  greatest  emphasis 
was  placed  on  examining  the  length  and  structure  of  the  transition  region  and 
the  parameters  which  influence  it.  The  second  phase  of  this  study  was  devoted 
to  examining  the  flow  over  "highly  ablated"  nose  shapes,  over  which  flow  sepa¬ 
ration  and  gross  flow  instabilities  occurred.  Here  we  studied  how  nose  tip 
configuration,  roughness  and  model  incidence,  as  well  as  the  Mach  number  and 
Reynolds  number  of  the  free  stream,  influenced  the  occurrence  and  structure  of 
large  scale  periodic  instabilities.  These  studies  demonstrated  the  strong 
effect  of  pressure  gradient  and  roughness  on  the  length  and  characteristics 
of  the  transition  region.  In  turn  these  parameters,  along  with  nose  tip  geo¬ 
metry,  strongly  influence  the  occurrence  of  flow  instabilities  over  highly- 
indented  nose  shapes. 

In  PART  II  of  this  investigation,  we  studied,  in  detail,  the 
generation  of  disturbances  in  the  stagnation  region  of  the  flow  as  a  particle, 
which  comes  from  the  model  surface,  interacts  with  the  bow  shock.  Here  we 
were  primarily  interested  in  gaining  a  knowledge  of  the  magnitude  and  mechanism 
of  heating  enhancement.  These  studies  which  were  conducted  at  Mach  6  and  13, 
demonstrated  that  there  was  little  heating  enhancement  as  the  particle  (between 
100  and  800  microns  in  size)  transversed  the  shock  layer;  however  heating  levels 
between  3  and  10  times  the  stagnation  point  value  can  be  generated  as  the  par¬ 
ticle  penetrates  the  bow  shock.  Four  distinctly  different  flow  regimes  were 
found  to  exist  for  different  penetration  distances  and  particle  trajectories. 

Of  these  the  most  dramatic  is  one  where  a  flow  instability  similar  to  that  en¬ 
countered  over  the  highly  indented  nose  shapes  is  observed. 
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EXPERIMENTAL  STUDIES  OF  TRANSITIONAL  BOUNDARY  LAYERS 
AND  BOUNDARY  LAYER  SEPARATION  WITH  GROSS  FLOW 
INSTABILITIES  OVER  ABLATED  NOSE  SHAPES 
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Section  1 
INTRODUCTION 


To  develop  a  reliable  and  accurate  re-entry  system  employing  slender 
high  /3  vehicles,  it  is  necessary  to  understand  and  formulate  methods  to  analyze 
the  changes  which  occur  to  the  ablative  nose  tip  and  the  associated  viscous  and 
inviscid  flow  field  in  which  it  is  enveloped  during  re-entry.  Small  changes  to 
the  shape  of  a  nose  tip  can  change  significantly  the  drag  and  moment  coef¬ 
ficients  of  the  vehicle,  particularly  if  such  shape  changes  are  asymmetric  in 

character.  While  measurements  from  ablation  tests  have  suggested  that  nose 

* 

shapes  develop  equilibrium  or  self-preserving  profiles  under  fixed  free  stream 
conditions,  complex  non- equilibrium  shapes  can  develop  as  the  vehicle  descends 
through  the  atmosphere,  and  its  boundary  layer  undergoes  a  transition  from 
laminar  to  turbulent  characteristics.  The  flow  fields  over  such  "non-equilibrium" 
shapes  are  complex,  involving  regions  of  shock  wave-boundary  layer  interaction, 
and  shock-shock  interaction.  The  large  localized  heating  rates  associated  with 
such  interactions  can  ultimately  create  nose  shapes  over  which  flow  instabilities 
resulting  from  the  alternate  formation  and  collapse  of  embedded  separated 
regions  can  occur.  The  d>Tiamic  forces  which  result  from  such  instabilities 
seriously  degrade  the  performance  of  the  vehicle.  Experimental  studies  have 
demonstrated  that  small  asymmetries  in  surface  roughness  or  ablated  shape  for 
these  configurations  can  result  in  flow  asymmetries  of  dramatic  proportions. 

A  small  angle  of  incidence  can  result  in  flow  asymmetries  of  comparable  magni¬ 
tudes  . 


The  movement  of  the  transition  point  from  the  conical  surface  onto 
the  face  of  a  nosetip  begins  a  process  which,  because  of  resultant  localized 
increase  in  heating,  may  create  an  indented  nose  shape.  The  formation  of  an  im¬ 
bedded  shock  over  this  configuration  results  in  a  significant  increase  in 
heating  at  the  base  of  the  shock  (resulting  from  shock  wave-boundary  layer  in¬ 
teraction)  which  drives  the  ablating  nose  into  a  shape  with  a  more  concave 
form.  On  these  shapes,  the  shock-shock  interaction  can  also  cause  intense 
heating  at  the  reattachment  point  of  the  shear  layer  which  is  generated  at  the 

* 

Discussed  in  Section  3. 
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intersection  between  the  bow  shock  and  embedded  shock.  In  many  cases  nose 
tips  do  not  ablate  symmetrically  and  the  asymmetric  flow  field  which  results 
can  impress  a  pressure  distribution  on  the  nose  which  can  cause  asymmetric 
forces  of  significant  proportions.  These  side  forces  are  of  greatest  magnitude 
when  a  separated  region  is  induced  on  one  side  of  the  nose  tip  while  the  flow 
remains  attached  on  the  other.  This  situation  can  also  arise  on  symmetric 
highly-ablated  nose  shapes  when  they  are  pitched  at  angle  of  attack. 

The  development,  through  ablation,  of  nose  shapes  over  which  gross 
aerodynamic  instabilities  occur,  has  been  demonstrated  in  ablation  tests  using 
camphor  models  and  inferred  from  downrange  observation  of  full  scale  flight 
missions.  Regions  of  gross  instability  on  nose  shapes  in  hypersonic  flow  were 
first  studied  in  detail  over  "spiked  bodies".  For  a  certain  class  of  configura¬ 
tions,  where  a  spike  pierces  the  bow  shock  and  yet  is  of  insufficient  length 
to  permit  an  equilibrium  separated  region  to  be  formed,  a  flow  oscillation  is 
observed  in  which  there  is  a  periodic  formation  and  collapse  of  the  separated 
region.  A  similar  oscillation  is  observed  over  bi-conic  configurations;  and 
studies  have  been  performed  to  define  the  geometries  over  which  these  oscilla¬ 
tions  occur.  Most  of  the  studies  which  were  conducted  in  hypersonic  flow  were 
at  low  Reynolds  numbers;  however  the  recent  high  Reynolds  number  studies  at 
Calspan  have  demonstrated  that  transition  to  turbulence  on  the  nose  of  the 
vehicle  can  significantly  inhibit  such  oscillations.  Nose  tip  roughness, 
which  is  an  inherent  feature  of  ablating  nose  shapes,  was  also  found  to  play 
an  important  role  in  the  oscillation  onset.  The  detailed  heat  transfer,  pres¬ 
sure  and  skin  friction  measurements  made  in  the  Calspan  study  identified  the 
importance  of  shock-boundary  layer  interaction  and  shock-shcok  interaction  as 
the  basic  mechanisms  which  cause  large  heating  rates  driving  an  indented  shape 
toward  shapes  over  which  flow  instabilities  occur.  Perhaps  the  most  important 
finding  from  these  studies  is  the  large  sensitivity  of  the  oscillating  flow 
fields  to  angle  of  attack.  The  large  asymmetric  pressure  forces  and  heat  transfer 
rates  generated  under  these  conditions  tend  to  drive  the  nose  region  into  an 
asymmetric  form  and  the  vehicle  to  increasing  angle  of  incidence.  The  oscil¬ 
lations  generated  under  these  conditions  differ  significantly  from  those  observed 
at  zero  incidence. 
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In  this  study  we  investigated  two  aspects  of  the  aerodynamics  of  the 
non- equilibrium  ablation  of  nose  tips.  In  the  first  of  these  studies  we  obtained 
detailed  heat  transfer  and  pressure  measurements  on  blunt  elliptic,  blunt 
biconic  and  triconic  models  at  Mach  number  from  8  to  13  for  Reynolds  numbers 
up  to  80  X  10^.  Here,  while  we  sought  information  on  the  laminar  heat  transfer 
distribution  over  these  nose  shapes,  we  were  primarily  interested  in  examining 
the  length  and  structure  of  the  transition  region  and  the  parameters  which 
influenced  it.  In  the  second  study  we  examined  the  flow  over  highly  ablated 
nose  shapes  over  which  flow  separation  and  gross  flow  instabilities  occurred. 

Here  our  studies  were  devoted  to  examining  how  nose  tip  configuration,  surface 
roughness  and  model  incidence,  as  well  as  Mach  number  and  Reynolds  number  of 
the  inviscid  flow  influenced  the  occurrence  and  structure  of  large-scale 
periodic  aerodynamic  instabilities  which  can  develop  over  these  shapes. 
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Section  2 

STUDIES  OF  TRANSITION  AND  RE-LAMINARIZATION  OVER  MILDLY  ABLATED  NOSE  SHAPES 
2.1  imODUCTION 

The  development  of  computer  codes  to  describe  accurately  the  shape  change 
of  a  nose-tip  of  a  slender  RV  as  it  descends  through  the  atmosphere  revolves 
about  an  accurate  description  of  the  boundary  layer  and  its  interaction  with  the 
ablation  material.  While  the  description  of  the  development  of  nose  tip 
shapes  in  laminar  flows  may  be  considered  to  be  relatively  accurate,  the 
complexities  in  the  fluid  mechanics  which  arise  when  transition  begins  on  the 
bodies  has  yet  to  be  fully  understood  or  successfully  described  in  an  analytical 
or  numerical  framework.  It  is  clear  that,  at  this  time,  in  order  to  assist  in 
understanding  the  origins  and  development  of  turbulence  over  complex  nose  shapes, 
we  require  an  insight  through  experimental  measurements.  The  problem  of  pre¬ 
dicting  the  position  of  boundary  layer  transition  on  ablating  nose  shapes  has 
yet  to  be  solved.  While  this  problem  is  difficult  enough  for  smooth  bodies,  on 
nose  tips  we  are  faced  with  understanding  the  influence  of  ablation  and  rough¬ 
ness.  While  recent  analytical  studies  have  used  linearized  stability  theory 

1  2  3 
(Ko)  and  turbulent  energy  methods  (Wilcox  and  Finson  )  in  attempts  to  describe 

this  process  with  some  success,  each  employs  some  degree  of  empiricism  in  their 

approaches.  In  correlating  measurements  of  transition  made  in  wind  tunnels,  we 

are  faced  with  the  problem  of  evaluating  the  influence  of  tunnel  noise  on  this 

process,  and  some  theories  (e.g.,  Wilcox's)  attempt  to  account  for  this. 

Aerotherm  have  suggested  from  their  studies  that  transition  can  be  correlated 

in  the  form  ’2'^*  Where  -t  is  the  roughness  height  and  0  the 

momentum  thickness  at  the  beginning  of  transition.  While  this  criteria  is 

simple,  it  has  the  merit  of  being  in  relatively  good  agreement  with  measurements 

made  by  Aerotherm  at  Mach  5.  Measurements  of  transition  over  sharp  flat  plates 

and  cones  in  the  Calspan  Shock  Tunnels  at  high  Mach  numbers  have  been  found  to 

be  in  remarkably  good  agreement  in  the  ballistic  range  and  correlations  of 

flight  data  as  shown  with  measurements  in  Figure  1. 

While  the  beginning  of  transition  may  be  influenced  by  tunnel  noise, 
the  characteristics  of  the  transitional  region  and  the  turbulent  boundary  layer 
which  developed  downstream  of  it  does  not  exhibit  this  sensitivity.  The  length 
and  characteristics  of  the  transition  region,  however,  do  depend  strongly  on  the 
pressure  and  entropy  gradients  in  the  inviscid  flow  and  the  roughness  of  the 


5 


TRANS 


surface.  For  constant  pressure,  high  speed  boundary  layers  in  the  absence  of 
entropy  graident,  the  length  of  the  transition  region  can  be  related  to  the  length 
of  the  laminar  run  as  shown  in  correlations  of  measurements  on  flat  plate  and 
cones  made  at  Calspan  presented  in  Figure  2,  However,  as  shown  in  Figure  3, 
the  length  of  the  transition  region  on  blunted  nose  tips  cannot  be  related 
simply  due  to  the  length  of  laminar  run.  Adding  roughness  not  only  increases 
the  level  of  heating,  but  also  significantly  modifies  the  length  and  character¬ 
istics  of  the  transition  region  as  shown  in  Figure  4.  The  photographs  obtainc-  ’ 
in  the  experimental  studies  on  bi-conic  configurations,  show  on  Figure  5 
indicate  a  wave- like  structure  in  the  boundary  layer  close  to  the  transition 
point.  This  structure  has  also  been  observed  in  experimental  studies  conducted 

4 

by  Demetriades.  The  records  of  the  heat  transfer  in  the  transition  region  on 
this  bi-conic  configuration  in  the  present  studies  show  distinctive  bursts  in 
heating  which  raise  the  local  heat  transfer  rate  to  close  to  the  turbulent 
value.  These  bursts  appear  to  be  convected  downstream  over  adjacent  gages. 

This  result  is  not  inconsistent  with  Emmons'  spot  theory  in  which  the  heat 
transfer  to  the  wall  is  related  to  an  intermittancy  parameter  and  through  the 
relationship 


Ch 


ic„ 


run 


+ 


ii-i) 


where  2^  is  the  fraction  of  time  the  boundary  layer  exhibits  turbulent  charac¬ 
teristics,  and  Cu  and  Cu  are  the  heat  transfer  coefficients  for  full 
turbulent  or  laminar  flow  respectively.  Clearly,  assumes  values  from  0  to 
1  through  the  transition  region  and  can  be  related  to  the  position  in  the 


transition  region  by  the  expression  V  = 


i-e 

1-e' 


as  shown  in  Figure  6. 


A  is  a  parameter  chosen  from  correlation  of  the  experimental  data.  While 
this  general  approach  has  been  used  with  some  success  to  describe  the 
distribution  of  heat  transfer  to  the  wall  in  the  transition  region,  it  is 
clear  it  must  be  modified,  based  on  experimental  measurements,  when  transition 
is  influenced  by  roughness,  vorticity,  entropy  and  pressure  gradients.  It  is 
of  interest  to  note  that  while  the  Emmons'  theory^  is  in  relatively  good  agree¬ 
ment  with  the  mean  values  of  heat  transfer  in  the  transition  region  at  high 
Mach  numbers,  the  fluctuation  measurements  of  heat  transfer  made  under  these 
conditions  do  not  display  the  intermittent  laminar  and  turbulence  character 
upon  which  this  theory  is  based. 
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Figure  4a  EFFECT  OF  ROUGHNESS  ON  TRANSITION  AND  RELAMINARISATION  OVER  THE 
BICONIC  CONFIGURATION  AT  MACH  13 
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Figure  4b  EFFECT  OF  ROUGHNESS  ON  TRANSITION  AND  RELAMINARkWION  OVER 
THE  BICONIC  CONFIGURATION  AT  MACH  11 
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Figure  5  SCHLIEREN  PHOTOGRAPH  SHOWING  BOUNDARY  LAYER  TRANSITION  OVER 
A  BLUNT  BICONIC  MODEL  AT  MACH  11 
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Figures  HEAT  TRANSFER  DISTRIBUTION  IN  THE  TRANSITION  REGION 
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One  char.icteristic  of  the  transition  region  which  is  important  in 
nose  tip  technology  is  the  "heat  transfer  overshoot"  which  occurs  close  to  the 
end  of  the  transition  process.  The  generation  of  turbulence  in  the  transition 
region  is  such  that  the  turbulent  scale  size  larger  than  that  for  a  corres¬ 
ponding  fully  turbulent  boundary  layer  is  generated,  and  this  results  in  in¬ 
creased  heating  levels.  This  "low  Reynolds"  number  turbulent  structure  is  in¬ 
fluenced  by  viscosity  and  eventually  dissipation  drives  the  turbulent  scale 
sizes  toward  the  local  equilibrium  values.  However,  this  relaxation  process 
can  occur  over  a  region  as  long  as  50  boundary  layers  thicknesses,  as  shown  by 
Holden 6  and  Bushnell.^  While  this  low  Reynolds  number  effect  causes  heating 
rates  of  only  125%  in  excess  of  the  equilibrium  turbulent  levels,  sophisti¬ 
cated  shape  change  codes  which  are  now  being  developed  should  be  formulated  in 
such  a  way  that  the  flow  downstream  of  transition  is  modelled  correctly.  Again 
experimental  measurements  are  needed  to  inspire  and  evaluate  such  models. 


2.2  EXPERIMENTAL  STUDIES  OF  LAMINAR,  TRANSITIONAL  AND  TURBULENT  FLOWS, 

INCLUDING  THE  EFFECTS  OF  ROUGHNESS,  ON  ABLATED  NOSE  SHAPES 

The  three  models  used  in  the  investigation  of  transition  on  ablated 
nose  shapes,  chosen  as  representative  of  those  that  might  be  found  in  flight, 
are  shown  in  Fig.  7.  These  models  were  highly  instrumented  with  heat  trans¬ 
fer  and  pressure  gages  to  provide  the  spatial  and  temporal  resolution  required 
to  define  accurately  the  characteristics  of  the  transition  and  relaminarisa- 
tion  region.  Figure  8  shows  the  model  with  the  elliptic  nose  shape  installed 
in  the  96"  Shock  Tunnel.  Schlieren  photographs  typical  of  those  obtained  in 
the  studies  are  shown  in  Figs.  9,  12  and  15.  The  distributions  of  heat  transfer 
over  the  biconic,  elliptic  and  tri-conic  configurations  at  Mach  numbers  from  8 
to  13  for  a  range  of  Reynolds  numbers  from  10°  to  80  x  10°  are  shown  in 
Figures  10,  13  and  16.  The  pressure  distributions  over  the  bi-conic  and 
elliptic  configurations  varied  little  with  Mach  number  and  Reynolds  number,  and 
typical  distributions  obtained  on  these  configurations  are  shown  in  Figs.  11 
and  14.  The  pressure  measurements  obtained  on  the  tri-conic  configuration  are 
shown  in  Fig.  17.  Here  the  variations  of  pressure  or  heating  distribution  with 
Mach  number  and  Reynolds  number  reflects  changes  in  the  structure  of  the  regions 
of  viscous  interaction  and  flow  separation  over  the  convex  surface  of  the  model. 
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Figure  7  MODELS  OF  THE  MILDLY  ABLATED  NOSE  SHAPES 
USED  IN  THE  EXPERIMENTAL  STUDIES 
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Table  1 

TEST  CONDITIONS  FOR  TESTS  OF  ABLATED 
NOSE  SHAPES  [SEE  FIGURE  7]  (CONTINUED) 
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Figure  8  ELLIPTIC  NOSE  TIP  MOUNTED  IN  96"  SHOCK  TUNNEL 
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Figure  9  A  SCHLIEREIM  PHOTOGRAPH  TYPICAL  OF  THOSE  OBTAINED  OVER 
THE  BICONIC  CONFIGURATIONS  (M  =  11) 
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Figure  10a  DISTRIBUTION  OF  HEAT  TRANSFER  OVER  THE 
(IVI-12.6,  Ren -2.1x10®) 


8  10  12 
BLUNT  BICONIC  CONFIGURATION 
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Figure  lOe  DISTRIBUTION  OF  HEAT  TRANSFER  OVER  BLUNT  BICONIC  CONFIGURATION 
(MACH  11  ReQ  =  8x  10°) 
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Figure  11a  TYPICAL  PRESSURE  DISTRIBUTION  ON  BICONIC  CONFIGURATION  AT  MACH  13 
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Figure  11b  TYPICAL  PRESSURE  DISTRIBUTION  ON  BICONIC  CONFIGURATION  AT  MACH  11 
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BICONIC 
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Figure  lie  TYPICAL  PRESSURE  DISTRIBUTION  ON  BICONIC  NOSE  AT  MACH  8 


Figure  12  SCHLIEREN  PHOTOGRAPH  OF  THE  FLOW  OVER  THE  ELLIPTIC  NOSE  TIP 
(M  =  13) 
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Figure  13b  MEASUREMENTS  OF  HEAT  TRANSFER  OF  THE  ELLIPTIC  NOSE  AT  MACH  11 
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Figure  14a 


A  TYPICAL  DISTRIBUTION  OF  PRESSURE  ON  THE 


ELLIPTIC  NOSE  AT  MACH  13 
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SOLID  SYM  120°  240°  RAYS 
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Figure  14b  A  TYPICAL  DISTRIBUTION  OF  PRESSURE  ON  THE  ELLIPTIC  NOSE  AT  MACH  7.8 
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Figure  15a 


A  SCHLIEREN  PHOTOGRAPH  OF  THE  FLOW  OVER  THE  TRI-CONIC  CONFIGURATION 
WITH  BOUNDARY  LAYER  SEPARATION  OVER  THE  SECOND  SURFACE  (M  =  13) 
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Figure  15b  SCHLIEREM  PHOTOGRAPH  SHOWING  BOUNDARY  LAYER  SEPARATION  FOR 
I  TRANSITION  FLOW  OVER  THE  TRI-CONIC  MODEL  (M  =  11) 
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Figure  15c 


SCHLIEREN  PHOTOGRAPH  FOR  TURBULENT  FLOW  OVER 
THE  TRI-CONIC  MODEL  (M  =  11 ) 
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Figure  16a  HEAT  TRANSFER  DISTRIBUTION  OVER  THE  TRICONIC  MODEL  AT  MACH  12.8 
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Figurs  17a  DISTRIBUTION  OF  PRESSURE  OVER  THE  BI  CONIC  CONFIGURATION  AT  MACH  12.8 
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Figure  17b  PRESSURE  DISTRIBUTION  OVER  THE  TRI-CONIC  CONFIGURATION  AT  MACH  11 


The  experimental  program  was  conducted  at  Mach  numbers  from  8  to 
13  for  a  range  of  Reynolds  numbers  such  that  laminar,  transitional  and  turbulent 
flows  were  obtained  over  the  model  nose  shapes.  Emphasis  was  placed  on  ob¬ 
taining  measurements  which  would  assist  in  attempts  to  evaluate  the  models, 
in  particular  those  of  the  transitional  and  relaminarization  regions,  which 
are  used  in  the  shape  change  codes.  We  performed  simple  calculations  to  de¬ 
termine  laminar  and  the  heat  transfer  distributions  over  the  biconic  and 

8  9 

elliptic  configurations  using  the  modified  Lees  equation  * 


= 


PcVc  Y' 


where  j  =  0  (2  dimensional),  j  =  1  (axi-symmetric) ;  /3  ^  J  ^  ^  ^  spheri¬ 

cal  nose  and  2  for  2:1  elliptic  nose.  In  this  expression  Pc  =  Pe/Po,  Uc  =  Ue/U®*, 
Y  =  rb/a  and  S  *  S/a,  where  a  is  a  reference  length  (a  =  Rjg  for  a  spherical 
nose  and  Rj^/4  for  the  2:1  elliptic  nose.  The  stagnation  pdnt  heat  transfer  was 
determined  using  the  Fay-Riddel 1  method.  For  turbulent  flow  we  use  the 
relationships  developed  by  Mirels  (REF  12) 


where 


t, 


REP 


/  O. 

\  D 


Pc  Vc 


:  0-1 


a/a 


02. 


)(^r 


and  AH  =  H  -  H  ,  k  =  Mirels'  roughness  coefficient  and  is  the  Prandtl 

0  w  R 

number.  The  three  comparisons  which  are  shown  in  Figs.  18,  19  and  20  demonstrate 
relatively  good  agreement  between  theory  and  experiment,  and  indicate  the  flow 
remains  completely  laminar  to  the  end  of  the  afterbody. 

The  characteristics  of  the  transition  region  over  the  bi -conic  con¬ 
figuration  were  explored  by  progressively  increasing  the  unit  Reynolds  number 
of  the  free  stream,  thereby  moving  transition  first  onto  the  conical  surface 
and  then  toward  the  hemispherical  nose  tip.  Typical  distributions  of  heating 
obtained  from  these  studies  are  shown  in  Fig.  21.  Here  we  see  than  the  length 
of  the  transition  region  is  not  simply  related  to  length  of  the  laminar  run 


46 


8 


RUN  21 


Figure  20  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT  FOR  LAMINAR  HEATING 
TO  THE  ELLIPTIC  NOSE  TIP 
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Figure  21  MEASUREMENTS  OF  TRANSITION  HEATING  ON  THE  BICONIC  CONFIGURATION 


ahead  of  transition  Oi  -  inomentum  thickness  at  the  beginning  of  the  transition 
process.  Our  measurements  indicated  that  the  transition  process  was  strongly 
influenced  by  the  local  pressure  gradient  on  the  inviscid  flow  and  the  entrain¬ 
ment  of  the  highly  non-honiogeneous,  rotational  inviscid  flow  -  the  entropy  layer 
originating  in  the  stagnation  region.  If  liowever  the  length  of  the  transition 
region  is  defined  from  experiment,  then  the  distribution  of  heat  tra’^sfer  through 
this  region  can  be  predicted  with  reasonable  accuracy  using  the  intermittency 
relationships  Cu  where  Y  was  dci'ined  earlier.  The 

comparison  between  the  simple  theory  and  experiment  sliown  in  Tig.  22  demonstrates 
that  the  mean  heat  transfer  rate  increases  monotonically  througli  the  transition 
region.  In  fact  the  temporal  variation  of  heating  at  a  specific  point  in  the 
transition  region  is  such  that  the  intermittency  model,  in  which  it  is  assumed 
that  the  surface  experiences  alternating  laminar  and  turbulent  conditions,  is 
supported  by  the  present  set  of  experimental  measurements. 

To  perform  accurate  calculations  of  the  boundary  layer  flow  over  the 
bi-conic  configuration  it  is  necessary  to  model  the  entrainment  of  the  entropy 
layer  into  the  boundary  layer,  i.e.  account  for  entropy  layer  swallowing.  A 
further  refinement  would  be  to  include  the  effects  of  vorticity  generated  in 
the  entropy  layer  into  the  model  of  the  boundary  layer  development.  To  assess 
the  effects  of  entropy  swallowing  it  is  desirable  to  perform  calculations 
to  determine  the  characteristics  of  the  inviscid  flow.  Such  calculations  have 
been  performed  by  Crowell^^  wlio  has  used  the  Godunov  method  to  determine  the 
inviscid  flow  over  the  elliptic  and  bi-conic  configurations.  The  Godunov  method 
was  used  by  Crowell  to  determine  the  shock  siiape  and  pressure  distribution  over 
the  elliptic  and  bi-conic  configurations.  The  comparisons  between  the  measured 
and  computed  shock  shape  shown  in  fig.  23  demonstrate  that  tlie  theoretical 
method  is  in  excellent  agreement  with  the  experimental  measurements. 

Model  roughness  had  a  significant  effect  on  the  magnitude  and  dis¬ 
tribution  of  heating  over  bi-conic  nose  shapes,  particularly  where  transition 
was  observed  on  the  "smooth"  configuration.  The  increase  in  heating  which 
resulted  from  introducing  0.005  roughness  on  the  45°  bi- conic  nose  shape  at 
Mach  13  is  shown  in  Fig.  24.  On  the  smooth  configuration  transition  began 
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Figure  24  HEAT  TRANSFER  MEASUREMENTS  DEMONSTRATING  THE  EFFECT  OF  5  MIL 
ROUGHNESS  ON  TRANSITION  AND  THE  CHARACTER  OF  THE  TRANSITIONAL 
AND  RELAMINARIZATION  REGIONS 
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on  the  conical  surface  approximately  2  inches  from  the  stagnation  point  and 
was  complete  by  the  cone-cone  junction.  The  output  from  the  heat  transfer 
gages  in  the  transition  region  exhibited  an  output  typical  of  the  intermittent 
character  of  the  turbulence  in  this  region,  where  the  timewise  variation  in 
heat  transfer  spanned  the  local  laminar  and  turbulent  values  (See  Fig.  24). 
Introducing  roughness  moved  transition  forward  to  a  point  just  downstream  of  the 
sphere-cone  junction  with  the  lieat  transfer  on  the  conical  surface  exhibiting  a 
fully  turbulent  character.  Our  measurements  demonstrated  that  transition 
and  the  heating  rates  in  the  turbulent  flow  were  symmetrically  positioned  on 
the  cone  surface.  It  is  clear  from  the  measurements  that  re-laminarization 
takes  place  at  the  cone-cone  junction  on  both  the  smooth  and  rough  body. 

However  while  re-transition  takes  place  on  the  rough  body,  the  flow  over  the 
afterbody  on  the  smooth  model  remains  laminar.  Clearly  re-laminarization  and 
the  effect  of  roughness  on  re-laminarization  are  phenomena  of  significant 
importance  in  nose-tip  technology. 

The  effect  of  roughness  on  the  transition  process  and  the  heat 
transfer  levels  downstream  of  transition  was  also  investigated  at  Mach  11  for 
free  stream  conditions  which  caused  the  beginning  of  transition  to  be  located 
close  to  the  midpoint  of  the  45°  bi-conic  surface,  and  at  the  junction  of 
this  surface  with  the  hemispherical  nose-cap.  These  measurements  are  shown 
in  Figures  25  and  26.  We  observe  that  roug  mess  has  moved  the  transition 
point  well  forward  for  the  tests  conducted  at  a  ARe^of  6  x  10^,  and  in  each 
case  has  inhibited  the  re-laminarization  process  which  occurred  at  the  sphere- 
cone  junction  on  the  smooth  bodies  .  While  relaminarization  occurred  on  the 
afterbody  downstream  of  the  cone-cone  junction  over  the  smooth  bodies  at  both 
Reynolds  numbers  the  flow  at  the  highest  Reynolds  number  on  the  rough  body 
remained  turbulent  through  this  expansion  region.  We  observe  from  the  measure¬ 
ments  made  on  both  the  smooth  and  rough  nose  tips  that  transition,  relarainariza- 
tion  and  the  distribution  of  heat  transfer  occurs  symmetrically  over  the  model 
surface. 


While  the  conventional  shape  change  cod:3  do  not  incorporate  features 
which  enable  the  position  of  re-laminarization  or  the  heat  transfer  distribu- 


tion  within  these  regions  to  be  calculated,  the  effects  of  roughness  has  been 

14 

incorporated  in  a  semi-empirical  technique.  Crowell,  using  the  Powars  roughness 

heating  formulation  together  with  shock  shape  and  pressure  distributions  from 

the  Godunov  code  obtained  the  calculations  shown  in  Fig.  27a  for  the  heat 

transfer  to  the  biconic  nose.  These  solutions  were  obtained  employing  normal 

shock  and  entropy  layer  swallowing  models  to  establish  the  local  properties 

outside  the  boundary  layer.  We  see  that  while  the  solution  employing  the  more 

attractive  entropy  layer  swallowing  technique  overpredicts  the  heating, 

that  based  on  the  normal  shock  approach  underestimates  the  heating  levels. 

It  should  also  be  noted  that  the  roughness  induced  augmentation  is  based  on 
14 

measurements  made  at  Mach  5,  and  the  good  agreement  between  the  measurements 
and  the  theory  in  the  absence  of  roughness  shown  in  Fig.  27b  suggests  that 
the  roughness  modelling  is  in  question.  Clearly  the  mechanics  of  roughness 
indue?  '  sition  and  heating  enhancement  and  wall  heating  in  regions  of 

transiti.  'd  re- laminari ration  remain  to  be  understood.  Experiments  are 

required  on  snooth  and  rough  nose-tips  over  a  range  of  bluntness  ratios  and 
free  stream  conditions  so  that  we  can  isolate  the  relative  effects  of  entropy 
layer  swallowing,  roughness  induced  turbulence  and  transition. 

The  relaminarization  of  a  turbulent  or  transitional  boundary  layer 
when  subjected  to  a  strong  streamwise  favorable  pressure  gradient  has  been 
observed  in  both  ground  and  flight  tests  and  is  clearly  an  important  phenomena. 
Relaminarization  not  only  causes  a  decrease  in  the  heating  levels,  but  creates 
a  situation  in  which  the  boundary  layer  can  be  more  easily  separated  by  an 
adverse  pressure  gradient.  This  latter  feature  is  of  critical  importance  in 
the  development  of  flows  over  highly  ablated  nose  shapes  where  boundary  layer 
separation  and  flow  instabilities  are  observed  to  occur.  The  large  decrease 
in  the  heating  levels  associated  with  re-laminarization  were  observed  in  the 
measurements  made  on  the  bi-conic  configurations  at  Calspan  and  are  shown  in 
Figures  25  and  26.  For  flows  without  roughness,  the  heat  transfer  gages  on 
the  conical  frustrum  always  give  mean  and  fluctuating  heat  transfer  values 
vdrich  were  in  close  agreement  to  those  for  fully  laminar  flow,  irrespective 
of  the  position  of  transition  on  the  fore-body.  Only  by  severely  tripping  the 
boundary  layer  with  roughness  is  one  able  to  preserve  their  turbulent  character- 
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istics  through  the  expansion  region,  as  shown  in  Figure  25.  While  some  criteria 
have  been  developed  to  suggest  the  conditions  under  which  boundary  layers  will 
re-laminarize  in  subsonic  flow,  no  such  criteria  exist  for  the  supersonic 
regime  with  roughness.  Relaminarization,  which  is  most  easily  understood  in 
terms  of  the  turbulent  energy  approach,  can  be  described  as  a  process  in  which 
strong  dissipation  of  the  turbulence  energy  takes  place  while  turbulence  pro¬ 
duction  is  reduced  because  the  favorable  pressure  gradient  decreases  the 
velocity  gradient  in  the  production  region  of  the  boundary  layer.  While  it 
may  be  possible  to  obtain  numerical  solutions  to  the  turbulent  energy  equations 
to  determine  the  development  of  turbulent  energy  in  highly  expanding  flows, 
experimental  measurements  are  required  to  provide  or  validate  the  numerical 
constants  used  in  such  an  approach. 


> 


61 


Section  3 


FLOW  SEPARATION  AND  THE  DEVELOPMENT  OF  INSTABILIEIES 
OVER  ABLATED  NOSE  SHAPES 


3. 1  INTRODUCTION 

Wind  tunnel  studies  using  ablating  models  have  demonstrated  that 
"slender"  concave  nose  shapes,  over  which  gross  flow  instabilities  are  observed, 
can  be  formed  under  certain  free  stream  conditions  and  nose-tip  properties. 

While  an  explanation  of  this  behavior  using  sophisticated  shape  change  codes 
awaits  the  development  of  accurate  models  of  regions  of  transition  and  relaminarisa- 
tion,  the  simple  equilibrium  solutions  obtained  by  Welsh^^  provide  an  important 
insight  into  the  development  of  ablated  nose  shapes.  Welsh's  analysis,  predicts 
the  occurrence  of  the  well  known  "blunt"  self-preserving  nose  shape  for  fully 
laminar  or  fully  turbulent  boundary  layers;  however  for  flows  where  transition 
occurs  ahead  of  the  sonic  point  on  the  nose  tip  (and  under  certain  fully  turbu¬ 
lent  conditions)  this  approach  predicts  the  occurrence  of  slender  nose  shapes 
similar  to  those  shown  in  Fig.  28.  This  result  has  been  verified  in  wind 
tunnel  tests  of  camphor  models  where  shapes  similar  to  those  shown  on  Figure  29 
can  be  developed  under  transitional  conditions.  The  existence  of  slender 
ablated  nose  shapes  have  been  demonstrated  in  many  experimental  studies  however 
in  general  these  shapes  exhibit  marked  asymmetric  characteristics .  Asymmetric 
nose  shapes  are  believed  to  result  from  asyinmetric  transition;  whatever  the 
cause  such  shapes  can  give  rise  to  large  de-stabilizing  moments  and  side  forces. 
Under  certain  conditions,  highly  idented  nose  shapes,  similar  to  that  shown 
in  Fig.  30  were  generated  on  camphor  models,  and  over  such  shapes,  the  viscous 
and  inviscid  flow  fields  interact  in  such  a  way  that  the  flow  became  oscillatory. 

In  studies  of  ablated  nose  shapes,  it  has  become  customary  to  describe 
the  nose  tip  geometry  in  terms  of  the  BI-CONIC  or  TRI-RADIUS  configuration. 
Schematic  representations  of  these  shapes  are  shown  in  Figure  31.  A  highly 
ablated  nose  shape  the  form  of  which  is  to  be  characterised  by  these  models 
is  shown  in  Fig.  32.  Bi-conic  configurations  with  ®  Qj  ~  ^  spiked 
bodies,  and  it  was  on  these  bodies  that  gross  instabilities  described  above 
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Figure  28  SLENDER  EQUILIBRIUM  SHAPE  DEVELOPED  UNDER  PURELY 
TURBULENT  BOUNDARY  LAYER  CONDITIONS  (FROM  REF.  10) 


0.49  tec  32.96  51.74  70.52  90.98  111.44 


Figure  29  VARIATION  OF  PROFILE  SHAPE  WITH  TIME  FOR 
AXISYMMETRIC  CAMPHOR  MODELS 
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17  1ft  20 

were  first  studied  in  hypersonic  flow  by  Maull,  Wood  and  Holden.  Maull 

and  Wood,  studied  the  flow  over  spiked  flat-ended  cylinders  and  cones  and  found 

(using  high  speed  schlieren  photography)  the  range  of  geometries  for  which 

gross  axial  pulsations  of  the  flow  are  observed  (see  Figure  33) .  This  flow 

* 

instability,  termed  the  "E"  oscillation,  occurs  on  a  blunt  body  when  the 
spike  pierces  the  bow  shock  and  yet  is  of  insufficient  length  to  permit  an 
equilibrium  separated  region  to  be  established.  While  Maull  and  Wood  used  only 
schlieren  photography  to  detect  flow  oscillations,  Holden  with  the  aid  of  high 
frequency  thin  film  gages  was  able  to  detect  two  additional  modes  as  well  as 
the  variation  of  heat  transfer  levels  which  occurred  during  these  oscillations. 
The  D  oscillation  arises  from  the  symmetrical  radial  enlargement  and  contrac¬ 
tion  of  the  separated  region  as  depicted  in  Figure  34.  The  D  oscillation 
has  a  frequency  slightly  larger  than  the  E  oscillation.  On  configurations 
which  were  close  to  the  E  oscillation  boundary,  Holden  observed  a  marked 
"lateral"  (L)  oscillation  which  appeared  insensitive  to  model  incidence.  It 
was  speculated  that  this  motion  was  a  circumferencial  oscillation  resulting 
from  an  inability  of  this  flow  to  generate  a  stable  reattachment  point. 

Again,  this  oscillation  had  a  frequency  larger  than  the  longitudinal  E 
instability.  Each  of  these  oscillations  were  accompanied  by  large  variations 
in  the  heating  rate  to  the  face  of  the  model.  A  common  feature  of  the  E 
oscillations  observed  in  these  tests,  in  which  only  laminar  flows  were  studied, 
was  that  the  Strouhal  number  (■^)  appeared  independent  of  Mach  number  and 
model  geometry. 

One  of  the  most  interesting  results  from  the  Pant  series  I  tests in 
which  flow  oscillations  were  observed  on  ablating  camphor  models  in  ground 
tests  was  that  these  instabilities  appeared  to  be  predominatly  of  the  D  and 
L  type,  although  the  E  oscillation  was  also  observed.  Seldom  was  a  completely 
symmetric  ablation  shape  observed;  in  general,  flow  instabilities  occurred  on 
one  side  of  the  model  and  enhanced  the  asymmetry.  An  ablated  nose  shape 
typical  of  one  where  oscillations  were  observed  is  shown  in  Figure  30.  It  is 
clear  that  this  shape  could  be  defined  in  terms  of  either  a  Bi-Conic  or  Tri- 
Radius  configuration.  The  former  representation,  however,  lends  itself  to  a 
* 

A  body  which  will  support  a  detached  shock  in  the  absence  of  a  spike. 
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LONGITUDINAL  OSCILLATION  IN  TYPE  E  OVER  SPIKED  BODIES  Figure  33 
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"D"  OSCILLATION  OVER  SPIKED  BODY  Figure  34 
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similar  treatment  to  that  used  to  define  the  oscillation  boundaries  for 
spiked-body  flows.  Figure  35  shows  the  oscillation  boundaries  determined  from 
the  spiked  body  studies  and  tests  on  Bi -Conic  configurations  performed  by 
Aerotherm^^.  The  results  from  the  latter  study  reinforce  the  conclusions  from 
the  studies  by  Wood  and  Holden  that  ©2  and  1/D  are  the  major  parameters  con¬ 
trolling  the  occurrence  of  the  E  oscillation.  The  oscillation  boundaries  de¬ 
termined  by  Aerotherm  for  turbulent  flow, fall  within  those  found  for  laminar 
flows  for  Wood  and  Holden.  Increasing  0^,  Rg/b  decreases  the  range  of 

the  and  1/D  for  which  flow  pulsations  are  observed.  The  effects  of  Mach 
number  and  roughness  were  not  included  on  the  Aerotherm  study.  Limited  variation 
in  Reynolds  number  was  included,  but  most  of  the  data  were  for  the  higher 
Reynolds  number,  turbulent  flow  condition.  While  the  studies  referenced  above 
reported  that  the  E  pulsations  for  highly  indented  nose  shapes  had  a  Strouhal 
frequency  of  close  to  0.2,  there  is  strong  evidence  that  instabilities  other 
than  this  oscillation;  namely,  the  D  and  L  modes  have  frequencies  which  are 
larger  and  influenced  by  the  model  geometry.  From  the  studies  performed  by 
Calspan  during  the  past  year,  it  is  clear  that  the  Mach  number  and  Reynolds 
number  and  roughness  are  important  parameters  which,  through  their  effect  on 
transition  and  boundary  layer  separation,  influence  the  occuirence  of  separation 
and  the  boundaries  between  the  D  and  E  oscillation. 

3.  2  STUDIES  OF  LAMINAR  FLOW  OVER  THE  WIDHOPF  1  AND  2  NOSE  SHAPES 

The  majority  of  the  studies  conducted  at  Calspan  were  performed  for 

conditions  such  that  the  boundary  layers  over  the  models  of  ablated  nose  shapes 

were  transitional  or  turbulent  in  character.  However,  we  studied  the  laminar 

20 

flow  over  two  ablated  nose  shapes  for  which  Widhopf  and  Victoria  obtained 
detailed  numerical  solutions  to  the  complete  flow  field  using  the  time- 
dependent  Navier-Stokes  equations.  These  nose  shapes  termed  the  1-10-10  and 
1-10-1  in  the  tri- radius  terminology  (or  tlie  Widhopf  1  ana  2)  are  shown  in 
Figure  36.  Detailed  heat  transfer  and  pressure  measurements  were  made  over 

both  of  these  configurations  at  Mach  10  for  a  Reynolds  number  (Re^^)  of  300. 

The  instrumentation  and  recording  equipment  were  designed  to  enable  frequencies 
up  to  1  Mz  to  be  recorded  with  little  attenuation.  Widhopf 's  calculations  for 
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STABILITY  BOUNDARY  FOR  E  OSCILLATION  ON  BICONIC  CONFIGURATIONS 


Table  2 

TEST  CONDITIONS  FOR  WIDHOPF  1  AND  2  CONFIGURATIONS 
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Rfc/FT.  -  3.708F+04  3.331E+-04  2.617E-+04  3.553E  +  04  3.4Z4E+04 

PITOT  PSIA  5.492E-01  5.506E-01  3.543E-01  5.590E-01  5.67bEHDl 

HIWI  FT  LBS/SLUG  3«186t+06  3.19bE*06  3,189E*-0t  3,183E+06  3.ia3E*-0fc 

TCrt)  °R  5,30bE+02  5.320E+02  5.310E^02  5.300E+02  5.300E*02 


the  1-10-10  configuration  showed  a  liigh  frequency  instability  occurred  in  the 
flow  field  close  to  the  point  of  surface  inflexion  on  this  body.  This  localized 
oscillation  was  jiredicted  to  have  a  major  frequency  of  approximately  600  KHz. 

Thus,  the  heat  transfer  gages  were  constructed  so  that  their  physical  dimensions 
were  smaller  than  the  wavelength  of  the  oscillation.  The  heat  transfer  distri¬ 
bution  to  this  body  is  shown  in  Figure  .T7  together  with  the  predictions  from 
the  Navier-Stokes  solutions  of  IVidhopf.  Although  the  measurements  are  in  good 
agreement  with  the  mean  values  from  these  calculations,  we  were  unable  to  de¬ 
tect  the  flow  instabilities  observed  by  Widhopf  in  his  numerical  solutions. 
Measurements  were  made  at  a  different  Mach  number  and  Reynolds  number  to 
determine  whether  an  oscillation  could  be  induced  at  slightly  different  flow 

conditions,  but  again  we  found  a  stable  flow.  However,  one  of  the  most  in¬ 
teresting  aspects  of  the  aerodynamics  of  this  configuration  is  the  large  heating 
rates  generated  at  the  inflexion  of  the  concave  surface  at  the  base  of  the 
embedded  shock  (see  Figure  38).  Studies  of  shock  wave-turbulent  boundary 
layer  interaction  on  concave  bodies  suggest  that  an  even  larger  increase  in 
heat  transfer  would  occur  on  the  1-10-10  configuration  if  the  boundary  layer 
were  transitional  or  turbulent.  An  ablating  body  upon  which  these  heating 
levels  were  impressed  would  tend  to  become  more  highly  indented  and  possibly 
develop  onto  the  shape  of  the  1-10-1  configuration. 

Heat  transfer,  pressure  and  schlieren  records  from  measurements  made 

on  the  1-10-1  configuration  at  Mach  10  demonstrate  that  the  flow  is  fully 

pulsating  in  the  E  oscillation  mode.  The  schlieren  photographs,  shown  in 

Figure  39,  indicate  the  oscillation  is  symmetric  with  a  frequency  which  gives 

fir 

a  Strouhal  number  (  )  ot  0.19.  The  heat  transfer  and  pressure  to  the 

model  (shown  in  Fig.  40)  undergoes  a  cyclic  variation  with  the  heat  transfer  vary¬ 
ing  between  a  value  which  is  approximately  equal  that  to  the  stagnation  point  on 
a  flat  faced  cylinder  with  the  same  body  diameter,  to  a  level  six  times  as  large. 
This  later  condition  is  attained  during  the  part  of  the  cycle  where  the  separated 
shear  layer  impinges  onto  the  almost  vertical  fact  'jf  the  afterbody.  A  small 
angle  of  attack  can  induce  a  second  asymmetric  uscill.\tional  mode  during  which 
a  circumferential  movement  of  the  separated  region  occurs.  This  mode,  which 
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Figure  37  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT  FOR  THE  WIDHOPF  1 
CONFIGURATION  (MO-10) 
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Figure  38  SCHLIEREN  PHOTOGRAPH  OF  FLOW  OVER  WIDHOPF  1  CONFIGURATION 


Figure  39a  SCHLIEREN  PHOTOGRAPH  OF  FLOW  OVER  WIDHOPF  2  CONFIGURATION 


Figure  39b  SCHLIEREN  PHOTOGRAPH  OF  FLOW  OVER  WIDHOPF  2  CONFIGURATION 


Figure  40b  HEAT  TRANSFER  AND  PRESSURE  MEASUREMENTS  ON  THE 
WIDHOPF  2  CONFIGURATION 
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will  be  discussed  in  greater  detail  as  it  occurs  both  on  the  "A"  shapes  and 
the  Pant  1  configurations,  induces  higher  harmonics  of  the  basic  E  oscillation, 
and  causes  asymmetric  loading  of  significant  proportions.  While  the  basic 
frequency  of  the  E  oscillational  mode  and  the  higher  harmonics  appear  to  be 
little  influenced  by  the  Mach  number  or  Reynolds  number  of  the  free  stream  or 
the  roughness  of  the  model  when  a  highly  ablated  nose  shape  is  tested,  for 
configurations  which  are  close  to  the  oscillation  boundaries  these  effects  do 
become  of  importance.  Such  is  the  case  for  a  series  of  "A"  nose  shapes  tested 
in  our  recent  experimental  program. 

3.3  EXPERIMENTAL  STUDIES  OF  THE  STABLE  AND  OSCILLATORY  FLOWS  OVER  THE 

"A"  SERIES  NOSE  SHAPES 


The  "A"  series  nose  shapes  were  developed  from  observations  made  on 
flight  and  ground  tests  to  establish  the  shape  of  ablated  nose  configurations  on 
which  pulsations  are  observed.  The  configuration  labeled  AO  in  Figure  41  was 
the  nose  shape  originally  chosen  by  AVCO  as  representative  of  a  configuration 
achieved  in  flight.  This  shape  was  tested  under  laminar  flow  conditions  and 
shown  to  exhibit  unsteady  flow  characteristics.  However,  configuration  A1  was, 
in  fact,  constructed  by  AVCO  to  be  tested  at  AEDC.  This  later  shape  was 
tested  at  Calspan  where  we  found  that  an  absence  of  E  pulsations  for  high 
Reynolds  number  conditions  at  Mach  11  and  13.  The  shapes  A2,  A3  and  A4  were 
selected  for  preliminary  studies  at  Calspan  to  provide  information  to  determine 
how  AVCO  could  make  simple  modifications  to  the  A1  configuration  to  obtain  an  E 
oscillation,  as  well  as  obtaining  basic  information  on  how  the  oscillation 
boundaries  are  modified  by  body  geometry,  Mach  number,  Reynolds  number  and 
roughness. 

The  characteristics  of  the  flow  field  over  the  A1  configuration  can  be 
seen  with  the  aid  of  the  high  speed  movie  sequence  shown  in  Figure  42.  A 
full  scale  E  oscillation  is  not  observed  on  this  configuration;  instead,  the 
D  oscillation  occurs.  Here,  there  is  a  periodic  growth  and  collapse  of  the 
separated  region,  however,  the  separation  point  remains  close  to  the  tip  of 
the  nose  rather  than  moving  from  the  tip  to  the  base  as  observed  for  the  E 
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Table  3 

TEST  CONDITIONS  FOR  THE  A  CONFIGURATIONS 
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Figure  42  THE  "D"  OSCILLATION  OVER  THE  A1 


oscillation.  The  mass  entrained  into  the  separated  region  is  partially  shed 
as  the  reattachment  point  in  the  separated  shear  layer  moves  from  the  shoulder 
of  the  nose  shape.  The  movement  of  the  shear  layer  toward  the  axis  accompanies 
the  subsequent  partial  collapse  of  the  recircular  region.  As  this  occurs,  the 
shear  layer  sweeps  across  the  face  of  the  nose  producing  heat  transfer  rates 
in  excess  of  five  times  that  to  the  stagnation  point  of  a  sphere  with  the 
same  base  diameter  as  the  nose  tip.  During  this  oscillation,  the  pressure  to 
the  nose  tip  undergoes  a  cyclic  variation  between  a  value  just  greater  than 
the  pitot  pressure  to  approximately  one  third  this  value,  while  the  nose  tip 
drag  varies  by  a  factor  of  approximately  two.  Further  measurements  with  the 
A1  configuration  suggested  that  the  D  oscillation  could  be  stabilized  with  in¬ 
creased  model  roughness,  Mach  number  and  Reynolds  number.  The  interesting 
feature  of  the  0  oscillation,  however,  is  that  a  Strouhal  number  of  0.30 
characterizes  the  frequency,  which  is  close  to  double  the  oscillation  frequency 
observed  for  the  basic  F  mode  oscillation.  In  fact,  our  measurements  suggest 
the  D  mode  may  form  the  first  harmonic  when  a  strong  E  oscillation  is  observed. 

Following  our  discovery  that  an  E  oscillation  could  not  be  induced 
over  the  A1  configuration,  studies  were  conducted  on  configurations  A2,  A3  and 
A4»which  were  shapes  derived  by  changing  the  tip  configuration  of  the  A1  shape. 
A2  has  the  same  bluntness  ratio  as  AO  but  a  slightly  larger  L/D.  We  found  that 
although  an  E  oscillation  was  generated  over  this  configuration  (see  Figure  43), 
this  oscillation  had  a  lower  non-dimensional  frequency  than  observed  earlier 
and  appeared  close  to  establishing  a  steady  separated  flow.  In  fact,  extending 
the  nose  to  a  sharp  point,  the  A3  configuration,  generated  a  steady  attached 
flow  (as  shown  in  Figure  44) ,  although  one  cycle  of  the  E  mode  was  observed 
just  before  flow  breakdown.  Here  we  believe  the  boundary  layer  ahead  of  the 
corner  radius  (R„)  was  fully  turbulent  and  for  this  reason  the  flow  did  not 

r 

separate  or  oscillate.  If  a  laminar  flow  were  to  be  established  over  the  tip 
of  tlie  A3  nose  shape,  we  believe  an  E  oscillation  would  result.  This  configu¬ 
ration  we  believe  would  be  sensitive  to  :  oughness  effects.  Finally,  our 
studies  of  the  flow  over  the  A4  configuration,  a  configuration  where  0^  was 
modified  to  achieve  approximately  the  same  bluntness  ratio  as  the  AO  configu¬ 
ration  was  found  to  be  just  unstable  in  the  E  mode,  as  shown  in  Figure  45. 
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Figure  43  MINI  "E"  OSCILLATION  OVER  THE  A2  CONFIGURATION 
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Thus,  nose  tip  configuration  Al,  A2,  A3  and  A4  have  geometries  which  induce 
flow  fields  which  lie  close  to  the  oscillation  boundaries  as  shown  in  Figure  25. 
However,  the  flow  over  the  AO  configuration  is  fully  pulsating  in  the  E  mode, 
and  this  is  the  configuration  over  which  detailed  heat  transfer  and  pressure 
measurements  were  obtained. 


3.4  DETAILED  MEASUREMENTS  ON  THE  AO  CONFIGURATION 


The  AO  configuration  was  fully  instrumented  with  high  frequency  re¬ 
sponse  pressure  and  heat  transfer  gages.  In  addition  to  placing  a  high  density 
of  instrumentation  along  the  0®  ray,  measurements  of  heat  transfer  and  pressure 
were  made  on  rays  inclined  at  120  and  240  degrees  to  this  principal  axis,  as 
shown  in  Figure  47.  In  this  way,  we  were  able  to  effectively  investigate  flow 
asymmetries  which  result  from  asymmetric  transition  or  angle  of  attack  effects. 

A  photographic  sequence  showing  the  flow  over  the  AO  configuration  at  zero 
angle  of  attack  is  shown  in  Figure  48.  Here  we  see  the  sequence  of  events  in 
which  the  bow  shocks  from  the  tip  and  the  afterbody  intersect  at  the  tip- 
afterbody  junction  producing  a  pressure  rise  which  the  boundary  layer  cannot 
sustain  without  separation  taking  place.  The  separation  point  moves  forward 
to  the  nose  tip  and,  because  it  is  impossible  to  establish  an  equilibrium 
separated  shear  layer  between  the  tip  and  the  shoulder  of  the  model,  the  sepa¬ 
ration  point  overshoots  the  nose.  The  resulting  flow  field  is  one  with  a  large  bow 
shock  standing  off  the  entire  configuration.  This  bow  shock  collapses  back  to 
the  model  surface  causing  the  entire  sequence  to  be  repeated.  Extremely  large 
heat  transfer  rates  are  generated  by  shock-boundary  layer  interaction  and  shock- 
shock  interaction  which  occur  at  the  junction  between  the  fore-  and  after-body. 
During  the  initial  formation  of  the  separated  region  at  a  time  when  a  strong 
oblique  shock  is  formed  at  the  forebody- afterbody  junction,  see  Figures  48 
and  49  severe  heating  is  observed  at  the  base  of  this  shock.  As  the  separated 
region  grows,  we  observe  the  formation  of  a  shear  layer  at  the  intersection 
between  the  fore-body  and  after-body  shocks.  The  heat  transfer  rates  which  are 
generated  in  the  region  where  this  shear  layer  attaches  to  the  surface  of  the 
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Figure  46  OSCILLATION  BOUNDARIES  ON  HIGHLY  INDENTED  NOSE  SHAPES 


H'-' 

(b)  A1  CONFIGURATION 


Figure  47  TWO  OF  THE  MODELS  FROM  THE  "A"  SERIES  TEST  PROGRAM 
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nose  tip  exceed  those  to  the  stagnation  region  of  a  sphere  with  the  same  diameter 
as  the  nose  and  after-body  by  factors  of  3  and  5  respectively.  The  typical 
outputs  of  heat  transfer  and  pressure  gages  positioned  in  the  junction  region 
arc  shown  in  Fig.  50.  We  see  the  heat  transfer  rate  varies  between  a  peak 
value  resulting  from  the  shock-boundary  layer  and  shock-shock  interaction  to  a 
minimum  value  corresponding  to  the  flow  configuration  where  the  bow  shock  is 
standing  off  the  nose  tip.  Further  distribution  of  heating  and  pressure  over  the 

AO  configuration  at  a  series  of  times  during  one  cycle  of  the  E  oscillation  is 
illustrated  in  Figs.  51  and  52.  It  is  clear  that  the  very  large  heating 
rates  at  the  center  section  of  the  model  would  lead  to  large  non-equilibrium 
distortions  in  the  shape  of  an  ablating  nose  tip.  The  structure  and  frequency 
of  the  oscillatory  flow  field  over  the  AO  configuration  showed  little  sensi¬ 
tivity  to  Reynolds  number  or  Mach  number  although  introducing  5  mil  roughness 
did  slightly  modify  its  amplitude.  However  the  supporting  of  this  model  at  an 
angle  of  attack  induced  a  flow  oscillation  which  was  asymmetric  in  character 
and  exhibited  two  distinct  modes  of  oscillation. 

The  asymmetric  flow  about  the  AO  configuration  at  an  angle  of  attack 
of  4.5°  is  shown  in  Fig.  53.  We  see  that  the  basic  mode  is  a  relatively  sym¬ 
metric  E  oscillation,  however,  during  the  period  when  the  flow  on  the  windward 
side  of  the  nose  tip  undergoes  one  cycle, the  flow  on  the  windward  side  under¬ 
goes  two  such  cycles.  Thus,  one  asymmetric  E  cycle  occurs  during  every  sym¬ 
metric  E  cycle.  Significant  asymmetric  pressure  and  heating  loads  are  created 
by  the  oscillation  as  indicated  by  the  distributions  shown  in  Fig.  54.  This 
asymmetry  is  a  maximum  at  the  time  when  the  flow  on  the  windward  side  is  at¬ 
tached  and  the  bow  shock  has  been  forced  to  move  just  beyond  the  tip  by  flow 
separation  on  the  leeward  side.  Thus, considerable  side  forces  are  generated 
which  in  addition  to  causing  a  pitching  momentum  at  the  nose  tip, exerts  severe 
bending  moments  which  could  possibly  lead  to  nose  tip  fracture. 


3.5  DETAILED  STUDIES  ON  THE  PANT  1  CONFIGURATION 


While  high  speed  movies  and  isolated  pressure  measurements  were  made 
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EFFECTS  OF  ROUGHNESS  ON  THE  PRESSURE  TRANSFER  MEASUREMENTS 
ON  THE  AO  CONFIGURATION  AT  MACH  13 


Figure  52b 
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Figure  53  ASYMMETRIC  OSCILLATION  ON  THE  AO  CONFIGURATION  AT 
3.50  ANGLE  OF  ATTACK 
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Figure  53  (Continued) 
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Figure  54a  EFFECTS  OF  ANGLE-OF  ATTACK  ON  THE  HEAT  TRANSFER  MEASUREMENTS 

ON  THE  AO  CONFIGURATION  AT  MACH  13  (Continued) 
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Figure  54b  EFFECTS  OF  ANGLE-OF-ATTACK  ON  THE  PRESSURE  TRANSFER  MEASUREMENTS 
ON  THE  AO  CONFIGURATION  AT  MACH  13 


during  the  PANT  test  series  conducted  by  Aerotherm  there  were  no  detailed 
measurements  of  the  heat  transfer  and  pressure  obtained.  Of  the  various  blunted 
bi-conic  configurations  studies  in  the  PANT  test  matrix,  the  PANT  1  configuration 
perhaps  best  represents  the  geometry  of  highly  ablated  nose  shapes  obtained  in 
ablation  studies  where  gross  instabilities  were  observed  in  the  inviscid  flow. 
Preliminary  studies  at  Calspan  had  also  revealed  that  under  some  circumstances 
significant  flow  asymmetries  can  develop  over  the  PANT  1  configuration  similar 
to  those  described  above,  for  the  AO  configuration  at  angle  of  attack.  To  de¬ 
tect  flow  asymmetries,  pressure  and  heat  transfer  instrumentation  was  placed  around 
the  model  circumference  as  well  as  along  one  major  instrumentation  ray.  A 
photograph  of  the  model  and  instrumentation  is  shown  in  Fig,  55.  The  studies 
over  this  configuration  were  conducted  at  Mach  11  and  13  over  a  range  of  Rey¬ 
nolds  numbers  so  that  laminar,  transitional  and  turbulent  flows  could  be  investi¬ 
gated.  A  sequence  of  photographs  taken  from  movies  of  the  flow  at  Mach  11  §  13 
are  shown  in  Figs.  56,  57  and  58.  We  observe  a  full  F.  oscillation  for  conditions 
where  we  believe  the  boundary  layer  is  transitional  over  the  body.  Adding 
roughness  did  not  markedly  affect  the  major  features  of  the  oscillation,  although 
some  degree  of  flow  asymmetry  is  induced  presumably  by  asymmetric  transition.  The 
major  fluid  dynamic  mechanisms  responsible  for  the  large  heating  rates  observed 
on  the  model  can  be  visualized  with  the  aid  of  Figs.  56-62,  the  shock-shock 
interaction  between  the  bow  shocks  from  the  forebody  and  the  afterbody  creates 
a  shear  layer  which  impinges  on  the  face  of  the  model.  Originating  close  to 
the  cone-afterbody  junction,  this  shear  layer  is  swept  outward  as  the  recirculation 
region  formed  over  the  cone  drives  the  cone  shock  away  from  the  surface.  Fi¬ 
nally  the  shear  layer  is  driven  to  the  afterbody  shoulder  and  the  air  trapped 
in  recirculation  region  formed  between  it  and  the  now  bulbous  cone  shock 
escapes  rapidly.  The  bow  shock  collapses  back  to  the  afterbody  surface  and 
the  entire  sequence  begins  again.  While  at  Mach  13  this  sequence  of  events 
is  not  significantly  modified  by  roughness,  at  Mach  11  introducing  5  mil  rough¬ 
ness  causes  the  separation  point  to  be  fixed  close  to  the  tip  of  the  forebody. 

The  flow  instability  which  is  observed  under  this  condition  resembled  a  strong 
"D"  type  oscillation  similar  but  more  violent  than  the  flow  oscillation  over 
the  A1  configuration.  Again,  we  observe  that  the  frequencies  associated  with 


109 


Figur«  55  PANT  1  CONFIGUr.ATION 
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Figure  56  FULL  E  OSCILLATION  OVER  SMOOTH  PANT  1  CONFIGURATION 
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Figure  58  FLOW  ASYMMETRIES  IN  THE  FLOW  ABOUT  ROUGH  PANT  SHAPE 
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these  modes  of  oscillation,  are  similar,  but  variations  of  pressure  and  heating 
rate  though  large  are  smaller  than  those  measured  when  a  full  1:  mode  is  observed. 


The  detailed  heat  transfer  distributions  on  the  Pant  I  show  that 
heat  transfer  rates  up  to  three  times  larger  than  the  stagnation  point  heating 
levels  to  the  nose  tip  of  the  model,  or  twelve  times  the  heating  rate  to  the 
stagnation  point  of  a  sphere  with  an  equal  base  diameter,  are  encountered  on 
the  second  conical  surface  of  the  model.  The  maximum  heating  rates  are  observed 
at  the  point  of  reattachment  of  the  shear  layer  formed  by  shock-shock  inter¬ 
action.  This  reattachment  point,  which  originates  just  above  the  cone-cone 
junction  is  swept  across  the  face  of  the  model  as  the  separated  region  grows, 
forcing  outward  the  position  of  the  shock-shock  interaction.  Our  interest  in 
the  L)  type  and  asymmetric  mode  instabilities  stem  from  the  fact  that  these 
phenomena  have  been  observed  in  ablation  tests  where  asymmetric  ablation  has 
induced  nose  shapes  which  are  at  most  equivalent  to  a  symmetric  nose  tip  placed 
at  angle  of  attack.  The  asymmetric  pressure  and  heat  transfer  and  pressure 
loading  can  play  a  key  role  in  the  survivability  and  accuracy  of  high  ^  re-entry 
vehicles. 


Section  4 
CONCLUSIONS 


A  detailed  study  has  been  conducted  of  the  aero-thermodynamics  of 
ablated  nose  shapes  at  Mach  numbers  from  8  to  13  and  Reynolds  numbers  up  to 
80  X  10^  based  on  nose  tip  diameter.  The  detailed  measurements  of  heat  transfer 
and  pressure  which  were  obtained  on  blunt  elliptic,  blunt  biconic  and  triconic 
nose  tip  configurations  demonstrated  the  importance  of  accurately  describing 
the  length  and  characteristics  of  regions  boundary  layer  transition  and  re- 
laminarization.  The  length  and  properties  of  the  transition  region  were  found 
to  be  strongly  influenced  by  the  local  pressure  gradient  and  the  structure  of  the 
entropy  layer.  However  while  the  length  of  this  region  (which  can  exceed  that 
of  the  preceding  laminar  lun)  cannot  be  predicted,  our  time-resolved  measurements 
support  the  intermittency  model  of  turbulence.  The  heat  transfer  measurements 
demonstrated  that  relaminarization  occurred  in  regions  of  strong  expansions  on 
all  the  nose  shapes  tested,  yet  this  phenomena  is  not  modelled  in  conventional 
shape-change  codes,  and  is  little  understood  as  a  basic  fluid  dynamic  effect. 

Model  roughness  was  found  to  have  a  major  effect  on  the  beginning  of 
transition  and  the  length  and  characteristics  of  the  transition  region.  Even 
when  the  beginning  and  length  of  the  transition  region  were  specified  from  ex- 
perminent,  theoretical  prediction  methods  based  on  the  Powers  roughness  model 
were  found  to  be  in  very  poor  agreement  with  measurements.  Roughness  does 
influence  relaminarisation;  however,  complete  relaminarisation  can  occur  in 
the  presence  of  roughness.  Regions  of  shock  wave-boundary  layer  interaction 
are  observed  on  mildly  indented  nose  shapes.  The  large  heat  transfer  rates 
generated  at  the  base  of  the  shock  formed  in  these  regions  of  viscous  inter¬ 
action  provide  the  mechanism  for  driving  ablating  bodies  into  more  highly  in¬ 
dented  nose  configurations. 

In  the  second  phase  of  this  investigation,  we  studied,  in  detail,  the 
structure  and  development  of  attached,  separated  and  unsteady  separated  flows 
about  highly  indented  nose  shapes.  While  emphasis  here  was  centered  on  the 
high  Reynolds  number  high  Mach  number  regime,  measurements  also  performed  on 
two  "highly  ablated"  nose  shapes  (Widhopf  1  and  2)  under  extremely  low  Reynolds 
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number  conditions,  =  300.  This  study  was  performed  explicitly  for  the 
purpose  of  comparing  high  frequency  heat  transfer  measurements  with  calcula¬ 
tions  made  with  the  Navier-Stokes  code  devised  by  Widhopf  and  Victoria.  While 
we  did  not  find  the  high  frequency  oscillation  predicted  by  Widhopf  and  Victoria 
over  the  Widhopf  1  configuration,  the  mean  heat  transfer  measurements  were  in 
excellent  agreement  with  theory.  The  measurements  on  the  Widhopf  2  configura¬ 
tion  exhibited  the  longitudinal  F.  mode  oscillation  with  a  characteristic 
Strouhal  number  of  0.17.  A  first  harmonic  of  this  oscillation,  an  asymmetric 
D  mode,  was  also  observed. 

Detailed  studies  were  made,  in  high  Mach  number  high  Reynolds  number 
flows,  of  the  occurrence  and  structure  of  flow  instabilities  occurring  on 
highly  indented  nose  tip  configuration  of  shapes  which  are  believed  to  be 
generated  under  some  flight  conditions.  In  the  first  set  of  studies,  we 
examined  the  model  configurations  and  free-stream  conditions  under  which  under¬ 
unsteady  flows  were  developed  over  highly  indented  nose  shape  under  turbulent 
flow  conditions  at  Mach  numbers  from  8  to  13.  The  oscillation  boundaries  were 
determined  for  smooth  highly-indented  bodies  and  expressed  in  terms  of  blunt-- 
ness  ratio  and  the  triconic  representation.  Measurements  on  selected  con¬ 
figurations  demonstrated  that  model  roughness  is  an  important  parameter  and 
should  be  explored  further.  E  and  D  oscillations  as  well  as  the  asymmetric 
L  mode  were  observed  in  these  studies.  The  detailed  heat  transfer  measurements 
over  the  AO  configuration  demonstrated  that  heat  transfer  rates  of  the  order 
ten  times  the  stagnation  point  value  to  a  sphere  of  an  equal  base  diameter  are 
generated  in  the  reattachment  region  on  highly  indented  nose  shapes.  Surface 
pressures  which  exceed  the  pitot  pressure  in  the  free  stream  are  also  observed 
in  this  region.  A  small  angle  of  attack  was  found  to  induce  strong  asymmetric, 
oscillations  which  could  induce  a  significant  destabilizing  moment.  The  sym¬ 
metric  E  mode  oscillation  was  found  to  have  a  Strouhal  number  between  0.15 
and  0.2  for  the  various  indented  nose  shapes  studies,  with  the  asymmetric 
oscillations  occurring  at  frequencies  which  were  harmonics  of  the  symmetric 
oscillations.  The  detailed  measurements  on  the  PANT  4  configuration  demon¬ 
strated  that  roughness  can  modify  the  E  mode  oscillation  and  induce  asymmetric 
oscillations  causing  asymmetric  forces  of  significant  proportions.  Again  ex¬ 
tremely  large  heat  transfer  rates  were  generated  in  the  reattachment  of  these 
flows. 
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1. 


INTRODUCTION 


The  phenomena  which  arise  when  the  nose  tips  of  hypersonic  re-entry 

vehicles  encounter  clouds  formed  from  ice  crystals  or  dust  particles  are 

known  to  be  of  key  importance  to  nose  tip  integrity  and  mission  accuracy.  Wliile 

nose  tip  erosion  resulting  from  particle  impact  has  been  a  primary  concern  of 

the  designer,  the  increase  in  nose  tip  heating,  relative  to  the  clean  air 

environment,  associated  with  particle  rebound,  or  the  ejection  of  material  from 

the  surface  of  the  nose  tip  can  play  a  critical  role  in  this  problem.  Recently, 

it  has  been  recognized  that  the  increased  heating  associated  with  simple  particle 

impact  and  associated  chemical  reactions  is  only  one  ingredient  of  the  augmenta- 

1  2 

tion  heating  problem.  Studies  in  dust  tunnels  at  Boeing  and  AEDC  have 
demonstrated  that  a  significant  increase  in  aerodynamic  heating  is  associated 
with  the  particle-flow  interaction  in  the  stagnation  region.  While  the  rebound¬ 
ing  particles  leave  a  wake  which  is  convected  to  the  surface,  changing  the  local 

2 

boundary  layer  properties,  studies  by  Wilkinson  et  al  at  TRW,  made  in  conjunction 
with  measurements  made  in  the  dust  tunnels  at  Boeing  and  AEDC,  have  identified 
bow  shock-particle  interaction  as  a  key  phenomenon.  Although  the  particles  must 
have  sufficient  momentum  to  reach  the  shock,  particle  size  does  not  appear  to 
strongly  influence  the  mechanism  of  heating  augmentation. 

Photographic  studies  made  in  the  wind  tunnels  at  Boeing  and  Calspan 
yielded  evidence  leading  to  the  identification  of  some  of  the  detailed  aero¬ 
dynamic  mechanisms  induced  by  particle-bow  shock  interaction.  Typical  photo¬ 
graphs  showing  a  single  particle-bow  shock  interaction  are  shown  in  Figures 
63  and  64.  Here  the  momentum  of  the  re-bounding  particle  is  such  that  only  a 
very  small  penetration  is  achieved.  Under  these  circumstances,  as  we  shall 
see  later,  a  ring  vortex  is  formed,  which  is  subsequently  convected  to  the 
surface  of  the  nose  tip  to  cause  increased  heating.  The  development  of  the 
vortex  ring  as  it  is  convected  toward  the  surface  was  demonstrated  first 
in  a  movie  sequence  (see  Figure  65)  taken  in  the  Boeing  tunnel.  The  non-planar 
behavior  of  the  vortex  ring  in  this  photographic  sequence  allows  us  to  observe 
its  expanding  toroidal  structure,  which  is  convected  to  the  surface  in  a 
fraction  of  a  millisecond.  While  particle- induced  vortex  generation  is  an 
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Figure  63  SHOCK-PARTICLE  INTERACTION  OVER  A  BLUNT  ELLIPSOID 
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Figure  64  SHOCK-PARTICLE  INTERACTION  OVER  A  BLUNTED  CONE 
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Figure  65  VORTEX  GENERATION  RESULTING  FROM  PARTICLE  REBOUND, 
A  SEQUENCE  TAKEN  IN  THE  BOEING  MACH  6  TUNNEL 


intriguing  phenomenon  which  can  cause  measurable  heating  augmentation,  it  is 
but  one  of  the  fluid  mechanical  mechanisms  which  result  in  particle- augmented 
heating.  IVhen  the  particle  has  sufficient  momentum  to  penetrate  through  the 
bow  shock,  a  shear  layer  or  jet  forms  at  the  intersection  of  the  particle  and 
bow  shocks  which  extends  to  the  model  surface,  presumably  causing  increased 
heating. 


Wliile  earlier  studies  in  dust  tunnels  demonstrated  significant 
particle-induced  heating  augmentation,  the  heat  transfer  and  pressure  instrumen¬ 
tation  lacked  the  frequency  response  to  follow  discrete  interactions.  Thus 
little  could  be  deduced  about  the  detailed  mechanisms  responsible  for  increased 
heating.  It  was  the  purpose  of  the  present  experimental  program  to  obtain 
time-resolved  measurements  of  heat  transfer  and  pressure  together  with  high  speed 
movie  photographs  to  establish  in  detail  a  quantitative  relationship  between 
particle-shock  interaction  and  surface  heating  augmentation.  In  the  following 
paper,  we  will  discuss  briefly  the  experimental  technique  that  was  employed, 
and  the  results  obtained,  in  the  studies  conducted  in  the  Cal  span  48  and  96- inch 
Shock  Tunnels.  Following  a  review  of  the  objectives  of  the  experimental  program, 
we  describe  the  design  of  the  experimental  program  and  the  specialized  instrumen¬ 
tation  developed  therein.  The  results  from  this  program  are  discussed  and  the 
detailed  heat  transfer  and  pressure  measurements  presented.  Here,  with  the  aid 
of  high  speed  schlieren  photography  we  demonstrate  that  a  number  of  dramatically 
different  flow  regimes  exist  depending  upon  the  mass  and  launch  velocity  of  the 
particle,  as  well  as  conditions  in  the  free  stream.  Heat  transfer  and  pressure 
measurements  for  these  interactions  are  presented.  The  conclusions  from  this 
study  are  given  together  with  recommendations  for  future  work. 

2.  EXPERIMENTAL  PROGRAM 

2.1  PROGRAM  OBJECTIVES 

The  objective  of  the  studies  conducted  at  Cal  span  was  to  delineate 
the  various  fluid  dynamic  mechanisms  responsible  for  heating  augmentation,  and 
determine  the  magnitude  of  this  augmentation  as  a  function  of  the  speed,  size 
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and  density  of  a  single  particle  ejected  from  the  surface,  the  geometry  of  the 
undisturbed  flow  and  the  Mach  number  and  Reynolds  number  of  the  free  stream. 

Here  heat  transfer  and  pressure  measurements  were  required  with  a  high  degree  of 
temporal  and  spatial  resolution  to  enable  specific  enhancement  mechanisms  to  be 
isolated.  We  first  wished  to  demonstrate  that  a  particle  fired  upstream  produces 
the  same  type  of  disturbances  observed  in  the  dust  tunnel  at  Boeing  where  the  dust 
was  added  to  the  airflow  close  to  the  throat  of  the  tunnel.  The  earlier  studies 
had  suggested  that  the  significant  heating  augmentation  occurred  only  when  the 
particle  possesses  sufficient  momentum  to  reach  the  bow  shock.  We  sought  to 
determine  the  magnitude  of  the  effect  on  nose  tip  heating  of  the  wake  shed  by 
the  particle  as  it  traversed  the  shock  layer.  By  launching  particles  of 
various  sizes,  densities  and  velocities,  we  wished  to  determine  whether  the 
resulting  particle-flow  interactions  were  influenced  by  particle  size.  For 
those  flows,  we  sought  to  establish  the  relationship  between  particle  tra¬ 
jectory,  'd  geometry  and  the  instantaneous  pressure  and  heat  distri¬ 

bution  over  1  .aodel.  Finally,  we  wished  to  determine  how  these  results 
were  influen  by  Mach  number  and  Reynolds  number. 

2.2  MODEL  AND  INSTRUMENTATION 

In  the  experimental  studies  performed  at  Cal  span,  we  chose  to  control 
the  particle  trajectory  by  launching  from  the  model  surface  rather  than  adding 
particles  to  the  free  stream  as  is  done  in  dust  tunnels.  While  this  approach 
gives  us  the  flexibility  and  control  required  to  perform  systematic  experi¬ 
mental  studies,  we  were  faced  with  the  problem  of  developing  the  sophisticated 
apparatus  required  to  launch  particles  between  100  and  800  microns  in  size  at 
velocities  between  30  and  300  ft/ sec.  A  major  task  in  this  research  program 
was,  therefore,  devoted  to  the  design  and  development  of  a  particle  launcher 
system  and  the  high  frequency  instrumentation  required  to  obtain  accurate 
time-resolved  measurements  of  heat  transfer  and  pressure  on  the  nose  shape. 

A  major  consideration  in  the  design  of  the  launcher  was  the  require¬ 
ment  that  the  launching  mechanism  caused  negligible  disturbance  to  the  air-flow 
in  the  stagnation  region,  or  create  mechanical  or  electrical  interference  which 
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would  disturb  the  output  from  the  delicate  high  frequency  instrumentation  on 
the  model  surface.  For  this  reason,  mechanisms  which  involved  a  direct  hy¬ 
draulic  launch  or  created  sizable  electrical  fields  were  discarded  in  our 
initial  design  studies.  After  numerous  approaches  had  been  explored,  we  chose 
a  technique  which  involved  impact  launching  of  the  particle.  A  diagram  of  the 
impact  particle  launcher  is  shown  in  Figure  66.  In  this  design,  we  employ 
two  pistons,  the  second  of  which  has  a  "cue"  mounted  in  the  center  of  its  face, 
which  is  placed  in  contact  with  the  particle  to  be  launched.  The  first  piston 
is  driven  down  the  launch  tube  when  the  driver  gas  (helium)  is  released 
behind  it  through  a  quick  acting  valve.  The  motion  of  this  piston  was  monitored 
with  photo  diodes  activated  by  LED's  mounted  along  the  internal  surface  of  the 
launch  tube.  The  second  piston  rests  in  contact  with  a  rubber  energy  absorption 
system  which  acts  to  restrict  its  travel,  when  struck  by  the  first  piston,  to 
within  0.50  inches.  The  particle  is  launched  when  the  elastic  compression  wave 
caused  by  piston-piston  impact  travels  down  and  is  reflected  from  the  face  of 
the  cue.  The  particle  travels  through  a  small  cavity  in  the  model,  the  geometry 
of  which  combined  with  the  expansion  wave  generated  as  the  second  piston 
rebounds,  produces  a  situation  in  which  the  pressure  disturbance  induced  by  the 
particle  as  it  exits  from  the  face  of  the  model  and  traverses  the  shock  layer 
is  effectively  the  only  disturbance  produced  in  the  launch  process.  Again, 
photo  diodes  and  LED's  are  used  in  a  measurement  scheme  to  determine  the  velocity 
of  the  second  piston.  A  specially  designed  inductive-capacitive  transducer 
placed  in  the  exit  plane  of  the  model  was  used  to  measure  directly  the  exit 
velocity  of  the  particles  when  the  launcher  was  calibrated. 

The  launcher  assembly  was  mounted  in  a  seismic  mass  which,  in  turn,  was 
supported  within  the  basic  model  by  a  soft  rubber  suspension  system.  In  this 
way,  the  instrumented  section  was  effectively  isolated  from  mechanical  distur¬ 
bances,  generated  when  the  launcher  was  activated.  The  instrumented  section  of 
the  model,  a  schematic  diagram  of  which  is  shown  in  Figure  67,  contained  80  heat 
transfer  and  28  high  frequency  pressure  gages.  Thin  film  resistance  thermometers 
were  used  to  obtain  detailed  temporal  measurements  of  surface  heat  transfer. 

These  gages,  which  were  fabricated  by  the  deposition  of  a  thin  platinum  film  on 
a  pyrex-substrate  provide  an  instantaneous  indication  of  the  surface  temperature. 

This  quantity  can  be  related  to  semi-infinite  slab  theory  to  the  instantaneous 

133 


134 


Figure  67  MODEL  ASSEMBLY  &  INSTRUMENTATION  NOSE  TIP  HEATING  (SAMSO/AFOSR) 


heat  transfer  rate.  Analog  networks  are  attached  to  each  gage  output  to  provide 
a  signal  which  is  directly  proportional  to  the  heat  transfer  rate.  Flush 
mounted  piezoelectric  pressure  transducers  with  a  frequency  response  up  to 
200  Kllz  were  used  to  obtain  the  temporal  variations  in  pressure  resulting 
from  the  interaction  between  particle  and  flow.  In  these  preliminary  studies 
a  flat-faced  cylindrical  shape  was  chosen  rather  than  a  hemisphere  because  it 
afforded  us  the  maximum  shock  layer  thickness  and  thereby  th';  maximum  resolution 
with  high  speed  cinematography.  The  basic  configuration  had  a  diameter  of 
6.0  inches;  however,  a  removable  ring  enabled  us  to  extend  the  diameter  to  9.0 
inches.  The  gages  were  positioned  so  that  a  high  density  of  instrumentation 
was  concentrated  in  one  quadrant  of  the  model,  while  the  balance  of  instrumen¬ 
tation  was  distributed  symmetrically  about  the  model  axis.  In  this  way,  we 
obtained  detailed  measurements  on  two  mutually  perpendicular  axes,  while  retain¬ 
ing  the  ability  to  track  asymmetric  disturbances  acrois  the  face  of  the  model. 
Both  steel  and  silica  spherical  balls,  with  dimensions  between  100  and  500 
microns  were  selected  to  be  launched  in  the  experimental  program.  The  size 
and  density  of  the  particle  and  the  velocity  of  the  launch  were  chosen  with 
the  aid  of  a  simple  numerical  code  so  that  the  particle  penetrated  specified 
distances  relative  to  the  position  of  the  shock  at  each  free  stream  condition. 
Mere  we  sought  to  understand  the  various  mechanisms  which  create  heating 
enhancement  for  the  different  flow  regimes  encountered,  as  well  as  determining 
whether  particle  size  ’s  an  important  factor  in  this  phenomena. 

2.3  EXPERIMENTAL  FACILITIES 

The  experimental  programs  were  conducted  in  Cal  span's  48-inch  and 
96-inch  Hypersonic  Shock  Tunnels  (Ref.  3  ).  The  operation  of  these  tunnels 
can  be  shown  simply  with  the  aid  of  the  wave  diagram  which  permits  the  high 
pressure  air  in  the  driver  section  to  expand  into  the  driven  section,  and  in 
so  doing,  generates  a  normal  shock  which  propagates  through  the  low  pressure 
air.  A  region  of  high  temperature,  high  pressure  air  is  produced  between  this 
normal  shock  front  and  the  gas  interface  between  the  driver  and  driven  gas, 
often  referred  to  as  the  contact  surface.  When  the  primary  or  incident  shock 
strikes  the  end  of  the  driven  section,  it  is  reflected  leaving  a  region  of 
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almost  stationary  high  pressure  heated  air.  This  air  is  then  expanded  through 
a  nozzle  to  the  desired  free  stream  conditions  in  the  test  section. 

The  duration  of  the  flow  in  the  test  section  is  controlled  by  tlie 
interactions  between  the  reflected  shock,  the  interface,  and  the  leading 
expansion  wave  generated  by  the  non-stationary  expansion  process  occurring 
in  the  driver  section.  At  Calspan,  we  normally  control  the  initial  conditions 
of  the  gases  in  the  driver  and  driven  sections  so  that  the  gas  interference 
becomes  transparent  to  the  reflected  shock,  as  shown  in  Figure  80a;  thus,  there 
are  no  waves  generated  by  interface-reflected  shock  interaction.  This  is  known 
as  operating  under  "Tailored-Interface"  conditions.  Under  this  condition,  the 
test  time  is  controlled  by  the  time  taken  for  the  driver-driven  interface  to 
reach  the  throat,  or  the  leading  expansion  wave  to  deplete  the  reservoir  of 
pressure  behind  the  reflected  shock;  the  flow  duration  is  said  to  be  either 
driver  gas  limited  or  expansion  limited,  respectively.  Figure  80b  shows  the 
flow  duration  in  the  test  section  as  a  function  of  the  Mach  number  of  the 
incident  shock.  Here  it  can  be  seen  that  for  operation  at  low  Ih,  running 
times  of  over  25  milliseconds  can  be  obtained  with  a  long  driver  section, 

3,  RESULTS  AND  DISCUSSION 

The  heat  enhancement  studies  at  Calspan  were  conducted  in  the  48"  and 
96"  shock  tunnels  at  Mach  numbers  from  6.5  to  16.3  and  Reynolds  numbers  from 
1  X  10^  to  30  X  10^/ft.  During  the  course  of  this  work,  22  runs  were  made;  a 
listing  of  these  runs  is  shown  in  Table  6 .  In  this  investigation,  we  made 
detailed  spatially  and  temporally  resolved  measurements  of  heat  transfer  and 
pressure  over  the  face  of  the  flat-ended  cylinder  shown  in  Figure  67.  The 
trajectory  of  the  particles  launched  from  the  model  were  varied  through  the 
choice  of  particle  size,  density  and  launch  velocity  to  examine  the  effect  of 
these  parameters  on  the  mechanism  and  magnitude  of  enhancement  heating.  The 
Mach  number  and  Reynolds  number  of  the  free  stream,  as  well  as  free  stream 
velocity  and  model  diameter,  were  varied  to  examine  how  transition  on  the 
model  and  in  the  free  shear  layer  influences  augmentation  heating.*  Here,  we 
also  sought  to  establish  whether  "Strouhal  number  scaling"  could  be  applied 
to  those  flows  in  which  gross  instabilities  were  observed  to  occur. 

k 

Test  conditions  given  in  Table  7 
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Table? 

TEST  CONDITIONS  FOR  HEATING  ENHANCEMENT  PROGRAM  (CONTINUED) 
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The  first  set  of  measurements  with  the  model  and  launcher  combination 
were  designed  to  calibrate  the  launcher  and  determine  the  extent  of  the  electri¬ 
cal  and  aerodynamic  effects  resulting  from  actuating  the  launch  mechanism. 

The  velocity  of  the  launched  particle  was  determined  as  a  function  of  the  driver 
pressure  and  piston  velocity  by  obtaining  the  lapse  time  as  the  particle 
traversed  between  two  inductive-capacitive  pickups  placed  coaxially  at  the  exit 
of  the  launcher.  With  the  sensitivity  of  the  recording  equipment  set  at  an 
order  of  magnitude  greater  than  required  at  the  lowest  heat  transfer  and  pressure 
levels  to  be  observed  we  found  there  were  no  disturbances  when  the  launcher 
was  actuated  firing  a  particle  in  a  vacuum.  Actuating  the  launcher  without  a 
particle  during  several  runs  demonstrated  that  there  was  no  measurable  difference 
in  the  heating  or  pressure  levels  before  and  after  activation  or  in  subsequent 
runs  at  the  same  conditions  where  the  launcher  was  not  operated. 

The  first  objective  of  the  tunnel  studies  was  to  define  the  types  of 

particle-flow  interaction  leading  to  heating  enhancement.  To  this  end,  a  series 

of  experiments  was  performed  adjusting  the  momentum  of  the  particle  so  that 

it  just  touched  the  bow  shock,  traveled  a  short  distance  beyond  it,  and  reached 

at  least  3  body  diameters  ahead  of  the  body.  Before  discussing  in  detail  the 

changes  in  flow  geometry  and  heat  transfer  and  pressure  distributions  which 

occurred  under  each  of  these  conditions,  it  is  important  to  comment  on  the 

measurements  while  the  particle  was  in  the  shock  layer.  For  the  range  of 

particle  sizes  and  launch  velocities  examined  (see  above)  there  was  no  measurable 

increase  in  the  heat  transfer  or  pressure  on  the  surface  of  the  model  as  the 

particle  traversed  the  shock  layer.  This  finding  is  in  agreement  with  the 

2 

results  of  Wilkinson,  et  al,  who  deduced  heating  augmentation  occurred  only 
when  the  rebounding  particles  possessed  sufficient  momentum  for  them  to  reach 
the  bow  shock. 

3.1  PARTICLE- INDUCED  VORTEX  HEATING  AUGMENTATION 

Heating  augmentation  was  first  observed  to  occur  for  trajectories  such 
that  the  particle  just  reached  the  bow  shock.  A  sequence  of  photographs  showing 
the  small  "dimple"  which  occurs  as  the  particle  reaches  the  bow  shock  and  the 
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subsequent  formation  of  a  toroidal  vortex  is  shown  in  Figure  68.  The  abrupt 
change  in  shock  curvature  at  the  junction  between  the  particle  and  body  shock 
creates  a  shear  layer  which  subsequently  rolls  up  into  a  vortex  ring  as  the 
particle  is  driven  back  into  the  shock  layer.  The  vortex  ring  is  then  convected 
toward  the  model  expanding  in  diameter  as  it  approaches  the  surface.  Correlating 
the  heat  transfer  time  histories  with  the  movie  sequence  reveals  that  increased 
heating  rates  are  observed  at  the  model  surface  over  a  period  of  just  over  1 
mil  1 isecond,  the  time  taken  for  the  vortex  to  traverse  the  shear  layer  and 
expand  across  the  model.  As  the  ring  passed  over  the  heat  transfer  gages, 
their  output  increased  from  1.5  to  3  times  their  undisturbed  values. 

3.2  HEATING  AUGMENTATION  RESULTING  FROM  SHOCK-SHOCK  INTERACTION  AND  FLOW 

SEPARATION 


Wlien  the  particle  has  sufficient  momentum  to  travel  ahead  of  the  ambient 
bow  shock  location,  very  dramatic  changes  occur  in  the  viscous  and  inviscid 
flow  fields  which  can  cause  heating  levels  that  exceed  the  abmient  by  factors 
of  2  to  4  close  to  the  model  axis  and  5  to  10  toward  the  periphery  of  the  disk. 
For  practical  purposes,  we  can  divide  the  phenomena  observed  into  three  classes: 
(1)  where  the  particle  travels  beyond  the  bow  shock  along  the  axis  of  symmetry 
but  reaches  less  than  0.7  body  diameters  from  the  model  surface,  (2)  where  the 
particle  travels  along  the  axis  of  symmetry  but  travels  beyond  0.7  body  diameter 
before  returning  toward  the  body  and  (3)  where  the  particle  travels  off  center 
setting  up  a  highly  asymmetric  flow  pattern. 

1.  Small  Upstream  Penetration  -  A  sequence  of  photographs  showing 
the  development  of  a  region  of  shock-shock  interaction  leading  to  heating 
augmentation  resulting  from  small  penetration  is  shown  in  Figure  69.  Here  the 
momentum  on  the  particle  was  just  sufficient  for  the  particle  to  reach  0.7  of 
the  body  diameter  ahead  of  the  model.  As  the  particle  moves  ahead  of  the 
bow  shock,  the  shock  wave  associated  with  the  particle  and  the  flow  behind  it 
interacts  with  the  original  bow  shock  inducing  a  shear  layer  at  their  point  of 
intersection.  This  shear  layer  sweeps  radially  outwards  across  the  face  of  the 
model  increasing  the  local  heat  transfer  at  its  base  by  as  much  as  a  factor 
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Figure  68 


heating  augmentation  associated  with  particle-induced 
RING  VORTEX  INTERACTION 
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of  ten.  In  this  sequence,  the  shear  layer  reaches  the  model  periphery  just  as 
the  particle  is  being  driven  back  to  the  model  surface,  and  the  sliear  layer 
then  collapses  inwards  in  a  quasi-steady  fashion,  disappearing  as  the  particle 
re-enters  the  shock  layer.  The  period  of  increased  heating  was  just  over  one 
millisecond,  the  time  taken  for  the  particle  to  exit  and  re-enter  the  shock 
layer.  Pressure  levels  just  above  the  pitot  level  were  observed  at  the  base 
of  the  shear  layer  while  those  within  the  annulus  corresponded  roughly  with  those 
achieved  by  compression  through  the  conical  shock  associated  with  the  particle. 

2.  Particle  Penetrations  Beyond  One  Body  Diameter  Resulting  in  Shock 
Layer  Oscillation  -  IVhen  the  particle  is  launched  with  sufficient  momentum  for  it 
to  reach  one  body  diameter  ahead  of  the  body  a  dramatic  change  occurs  in  both  the 
character  of  the  inviscid  flow  and  the  pressure  and  heat  transfer  on  the  model 
surface.  Here  the  entire  flow  field  begins  pulsating  as  demonstrated  in  the 
high  speed  movie  sequences  shown  in  Figures  70  and  71.  These  sequences,  which 
were  taken  at  Mach  6.5  and  13  respectively  at  similar  Reynolds  number  conditions 
demonstrate  considerable  similarity  in  mechanics  of  oscillation.  This  oscilla- 

*  4  r 

tion  IS  very  similar  to  pulsating  flows  observed  over  spikeu  bodies  or  highly 
indented  nose  shapes  as  demonstrated  by  the  movie  sequence  of  these  flows  shown 
in  Figure  48.  The  occurrence  of  these  massive  pulsations  appears  independent 
of  particle  size  and  once  the  particle  has  reached  a  body  diameter  ahead  of 
the  surface,  the  axial  motion  of  the  particle  is  strongly  influenced  by  the 
flow  oscillations.  The  non-dimensional  frequency  was  found  to  lie 

between  0.17  and  0.19  and  appeared  relatively  independent  of  Mach  number, 
particle  velocity  size  or  penetration,  and  model  size.  A  Strouhal  number  of 
between  0.17  and  0.20  typifies  the  frequencies  encountered  over  spiked  and 
highly  indented  bodies,  and  re-enforces  the  concept  that  the  basic  mechanism  is 
associated  with  a  simple  inviscid  filling  and  spilling  mechanism.  The  oscilla¬ 
tion  is  initiated  when  the  shear  layer  or  jet  formed  by  particle/ spike  shock-bow 
shock  interaction  re-attaches  to  the  body  surface  trapping  a  conical  region  of 
gas.  The  mass  which  is  being  continuously  added  to  this  region  through  the 
conical  shock  drives  the  local  stagnation  point  at  the  base  of  the  shear  layer 
toward  the  periphery  o‘f  the  model.  This  region  collapses  as  the  entrapped  gas 
escapes  when  the  shear  layer  moves  off  the  body  and  the  sequence  begins  again 
with  bow  shock  re-establishing  and  the  particle  shock-bow  shock  interaction 

f 

Termed  the  E  oscillation  in  this  study. 
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Figure  71 


(continued) 
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Figure  71  (continued) 


reforming  the  shear  layer.  Perhaps  the  most  intriguing  questions  is  why  two 
distinct  shock  structures,  one  around  the  minute  particle  and  a  second  around 
the  model  are  not  formed  at  this  juncture.  Both  the  pressure  and  heat  transfer 
records  exhibited  large  excursions  from  the  ambient  conditions.  A  typical  record 
■rom  a  thin  film  heat  transfer  gage  on  the  face  of  the  model  is  shown  in  Figure 
72  .  After  the  tunnel  starting  process  the  gage  is  exposed  to  a  period  of  constant 
heating  equal  to  that  observed  for  the  unperturbed  flow.  The  heating  rate 
increases  sharply  as  the  particle  interacts  with  the  shock  and  the  shear  layer 
sweeps  across  the  face  of  the  model.  A  decrease  in  heating  occurs  as  a 
bulbous  shock  layer  is  formed  and  the  flow  "collapses"  toward  the  body.  This 
pattern  is  repeated  in  a  flow  oscillation  of  surprising  regularity  and  persistence. 
While  the  hea*  ‘  ’  i-ates  across  the  model  vary  with  time,  the  magnitude  of  the 
maximum  vaTi.  ot  appear  sensitive  to  the  exact  particle  trajectory  as  can 

be  seen  from  ti  isurements  made  in  two  runs  (shown  in  Figure  73)  at  identical 
free  stream  conditions  but  with  different  particle  trajectories.  This 
figure  illustrates  that  the  edge  of  the  cylinder  experiences  the  largest  heat 
transfer  rates,  which  can  exceed  the  ambient  heating  levels  by  as  much  as  a 
factor  of  five.  The  minimum  heating  rates  fall  close  to  the  ambient  heating 
level  in  the  absence  of  a  particle  and  correspond  to  the  flow  condition  where 
the  bow  shock  is  collapsing  back  to  the  body.  A  typical  variation  of  the 
fluctuating  surface  pressure  close  to  the  center  of  the  model  for  this  condi¬ 
tion  is  shown  in  Figure  74.  The  maximum  levels  are  roughly  equal  to  the  pitot 
pressure  while  minimum  corresponds  to  the  plateau  pressure  for  conical  separated 
flows  over  spiked  bodies,  supporting  the  postulation  that  a  recirculation  region 
is  formed  over  the  model  during  part  of  each  cycle. 

The  basic  features  of  particle-induced  oscillations  at  Mach  12 
described  above  were  similar  to  those  observed  at  Mach  6.  A  photographic  se¬ 
quence  showing  the  E  oscillation  at  this  Mach  number  is  shown  in  Figure  48.  Tlie 
heat  transfer  distribution  for  this  condition  (shown  in  Figure  75J  again  indi¬ 
cates  that  heating  is  a  maximum  toward  the  edge  of  the  cylinder  with  an  augmenta¬ 
tion  factor  of  over  5;  however,  the  augmentation  factors  for  the  inner  gages  of 
approximately  3  are  slightly  larger  than  experienced  on  equivalent  gages  at  the 
higher  Mach  number.  Increasing  the  Reynolds  number  increased  the  magnitude  of 
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Figure  72  TYPICAL  HEAT  TRANSFER  RECORD  FROM  THIN  FILM  GAGES 
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Figure  73  HEAT  TRANSFER  DISTRIBUTION  TO  THE  FACE  OF  THE  DISC  FOR  A 
FULLY  OSCILLATING  FLOW  OVER  THE  MODEL 
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Figure  74  TYPICAL  PRESSURE  OSCILLOGRAM,  RUN  3 
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heating  to  the  center  section  of  the  cylinder  (see  Figure  76)  suggesting  that 
the  flow  is  transitional  or  turbulent.  Shown  on  Figures  75  and  76  are  measure¬ 
ments  of  the  maximum  heating  levels  observed  at  the  beginning  of  the  run  when 
minute  particles  of  the  mylar  used  to  seal  the  throat  of  the  nozzle  were 
convected  downstream  before  bouncing  from  the  face  of  the  model  and  interacting 
with  the  bow  shock.  These  shock-shock  interactions  produced  heating  rates 
which  were  comparable  with  those  induced  in  the  larger  scale  oscillation 
principally  because  the  same  mechanism  was  responsible  for  the  augmentation. 

3.  Augmentation  for  Off-Axis  Trajectory  -  Some  of  the  largest 
augmentation  factors  measured  in  the  present  studies  were  observed  for  trajecto¬ 
ries  where  the  particle  moved  off  axis  as  shown  in  the  photographic  sequence 
in  Figure  77.  Here,  because  of  the  asymmetric  nature  of  the  motion,  the 
mechanism  which  caused  flow  oscillation  for  the  degree  of  penetration  achieved 
by  this  particle  was  not  activated.  Instead  the  gas  which  normally  would  be 
trapped  within  a  recirculation  region  escaped  to  one  side  of  the  model  and  a 
non-oscillatory  interaction  region  swept  across  the  model.  The  maximum 
heating  levels  observed  on  the  face  of  the  model  are  shown  in  Figure  78.  These 
levels  were  observed  for  a  period  of  just  under  1  millisecond  and  the  increased 
heating  results  from  the  same  type  of  shock-shock  interaction  observed  for  the 
oscillatory  flows.  Finally,  the  photographic  sequence  shown  in  Figure  79 
provides  an  example  where  an  oscillation  is  first  observed  when  the  particle 
initially  travels  along  the  axis  of  symmetry;  however,  because  it  moves  off 
axis,  an  asymmetric  flow  field  results  in  the  spillage  of  air  from  one  side 
of  the  model  and  the  stabilization  of  the  oscillation.  The  particle  is  swept 
downstream  and  an  unperturbed  flow  re-established  over  the  model.  Augmentation 
factors  between  3  and  7  times  the  unperturbed  heating  levels  are  observed  for 
this  flow. 
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non-oscillatory  flow  resulting  from  off-axis  particle 
trajectory 
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Figure  77  (continued) 
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Figure  78  HEAT  TRANSFER  RATE  DISTRIBUTION,  RUN  20 
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3.3 


1-LUll)  DYNAMIC  MliCllANlSMS  IN  AUGMENTATION  lEATlNC. 


Studying  the  interaction  between  a  single  particle  and  the  shock  layer 
allows  us  to  focus  more  accurately  on  the  underlying  mechanisms  responsible  of 
heating  augmentation.  The  movement  of  a  single  small  particle  across  the  shock 
layer  can  be  predicted  with  reasonable  accuracy  using  simple  Newtonian  concepts 
employing  a  drag  coefficient  determined  for  free-molecular  flows.  In  our  studies 
with  steel  and  silica  particles  between  100  and  800  microns  in  size  we  were  unable  to 
detect  heating  augmentation  associated  with  interaction  between  the  wake  of  the 
particle  (in  tlie  shock  layer}  and  boundary  layer  over  the  model.  One  might  antici¬ 
pate  that  the  particle  wake  v^^ould  increase  the  vorticity  at  the  edge  of  the  bound¬ 
ary  layer  increasing  tlie  local  heating  and/or  "tripping"  the  laminar  boundary 
layer  in  the  stagnation  region.  If  such  effects  were  present,  they  did  not 
create  a  measurable  increase  in  heating  for  our  system  and  conditions.  The 
first  measurable  lieating  augmentation  occurs  when  the  toroidal  vortex,  formed 
when  the  particle  causes  a  minor  perturbation  of  the  bow  shock,  is  convected 
over  the  model  boundary  layer.  The  unsteady  changes  in  the  shock  curvature 
in  the  region  of  shock-shock  interaction  are  believed  responsible  for  the 
formation  and  subsequent  "roll-uji"  of  an  annular  shear  layer.  The  ring  vortex 
so  formed  expands  as  it  is  convected  over  the  body,  and  the  increased  heating 
can  be  explained  in  terms  of  increased  vorticity  of  the  edge  of  the  boundary 
layer.  If  this  vorticity  acts  to  "trip"  the  originally  laminar  boundary 
layer  then  the  heating  augmentation  can  be  correspondingly  larger.  When  there 
is  a  sustained  movement  of  the  particle  ahead  of  the  bow  shock  then  the  interac¬ 
tion  between  the  shock  surrounding  the  flow  field  associated  with  the  particle, 
and  the  how  shock  enveloping  the  model,  creates  a  shear  layer  or  jet  which 

attaches  to  the  face  of  the  model.  The  heating  rates  associated  with  such 

4  6 

interactions  have  been  explored  by  Holden  in  spiked  body  flows  and  Edney, 
and  Keyes  and  Mains  in  regions  of  shock  interference.  It  seems  clear 
that  major  augmentation  occurs  when  the  "Edney  Type  IV"  or  "jet"  interaction 
is  established  over  the  model.  Until  the  shear  layer  is  swept  to  the  edge  of 
the  model  for  the  first  time,  the  movement  of  the  particle  can  be  predicted  by 
simple  drag  models;  however,  when  the  flow  begins  to  oscillate,  the  alternate 
formation  and  collapse  of  the  recirculation  region  ahead  of  the  body  imposes  a 
system  of  forces  on  the  particle  which  makes  its  motion  extremely  complex.  The 
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particle  travels  far  beyond  the  point  predicted  by  applying  the  simple  laws  of 

motion  to  an  isolated  particle  in  the  free  stream.  As  discussed  earlier,  the 

mechanism  of  oscillation  is  controlled  primarily  by  simple  inviscid  flow 

r  r* 

considerations  and  the  non-dimensional  frequency  remains  remarkably 

uniform  at  0.16  and  0.19,  figures  similar  to  those  obtained  for  the  E  oscillation 
over  spiked  bodies.  Again,  the  maximum  heating  is  found  at  the  base  of  the 
re-attaching  jet  of  shear  layer  v.here  factors  of  3  to  6  times  the  ambient  un¬ 
perturbed  levels  are  observed,  i’he  minimum  heating  values  were  just  less  than 
those  measured  in  unperturbed  flow.  Pressures  on  the  flat  face  were  found  to 
vary  between  values  just  greater  than  the  pitot  pressure  to  minimums  roughly 
equal  to  the  plateau  pressure  corresponding  to  separated  flow  ever  the  face  of 
the  model. 

If  the  particle  travels  off  axis  permitting  an  asymmetric  interaction 
region  and  asymmetric  spillage,  then  a  non-oscillatory  interaction  region  can 
be  observed.  While  the  period  of  augmented  heating  is  less  for  these  cases, 
augmentation  factors  are  observed.  Again,  the  major  mechanism  is  interference 
heating  resulting  from  shock-shock  interaction. 

4.  CONCLUSIONS  AND  RECOMMENDATIONS 

An  experimental  study  has  been  presented  which  demonstrates  that  a 
single  minute  particle,  launched  through  the  bow  shock  of  a  flat-ended  cylinder, 
can  induce  severe  distortions  in  the  bow  shock  shape,  augmenting  the  heating 
rates  to  the  model  surface  by  as  much  as  a  factor  of  8.  This  investigation 
identified  four  distinctly  different  classes  of  flow  phenomena  resulting  from 
particle-shock  wave  interaction  together  with  temporally  and  spatially  resolved 
heat  transfer  and  pressure  measurements  associated  with  them.  The  most  dramatic 
of  these  results  when  the  particle  penetrates  more  than  a  body  diamel  jr  beyond 
the  bow  shock  and  massive  flow  instabilities  similar  to  those  encountered  over 
spiked  bodies  are  observed.  These  studies  have  identified  the  principle  mechanism 
of  enhancement  heating  as  one  associated  with  the  reattachment  of  a  jet/ shear 
layer  generated  by  particle  shock-bow  shock  interaction  after  the  particle  has 
penetrated  the  bow  shock.  Future  studies  should  investigate  how  nose  tip 
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configuration  influences  the  large  heating  augmentation  and  flow  instabilities 
observed  in  the  present  studies.  Measurements  should  be  made  to  determine 
particle-induced  augmentation  over  rough  models  under  fully  turbulent  boundary 
layers.  The  studies  should  be  extended  to  include  multiple  particle  interactions 
and  off-axis  launch.  Skin  friction  measurements  in  addition  to  the  measurements 
of  heat  transfer  and  pressure  should  be  made  to  better  define  the  flow  adjacent 
to  the  model  surface. 
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Appendix  A 

EXPERIMENTAL  FACILITIES  AND  MEASUREMENT  TECHNIQUES 

A.l  EXPERIMENTAL  FACILITIES 

The  experimental  programs  were  conducted  in  Calspan's  48-inch  and 
96-inch  Hypersonic  Shock  tunnels  (Ref.  14).  The  operation  of  these  tunnels 
can  be  shown  simply  with  the  aid  of  the  wave  diagram  shown  in  Figure  80a.  The 
tunnel  is  started  by  rupturing  a  double  diaphragm  which  permits  the  high 
pressure  air  in  the  driver  section  to  expand  into  the  driven  section,  and  in 
so  doing  generates  a  normal  shock  which  propagates  through  the  low  pressure 
air.  A  region  of  high  temperature,  high  pressure  air  is  produced  between 
this  normal  shock  front  and  the  gas  interface  between  the  driver  and  driven 
gas,  often  referred  to  as  the  contact  surface.  When  the  primary  or  incident 
shock  strikes  the  end  of  the  driven  section,  it  is  reflected  leaving  a  region 
of  almost  stationary  high  pressure  heated  air.  This  air  is  then  expanded 
through  a  nozzle  to  the  desired  free  stream  conditions  in  the  test  section. 

The  duration  of  the  flow  in  the  test  section  is  controlled  by  the 
interactions  between  the  reflected  shock,  the  interface,  and  the  leading 
expansion  wave  generated  by  the  non-stationary  expansion  process  occurring 
in  the  driver  section.  At  Calspan  we  normally  control  the  initial  conditions 
of  the  gases  in  the  driver  and  driven  sections  so  that  the  gas  interface  be¬ 
comes  transparent  to  the  reflected  shock,  as  shown  in  Figure  79  ;  thus,  there 
are  no  waves  generated  by  interface-reflected  shock  interaction.  This  is  known 
as  operating  under  "Tailored-Interface"  conditions.  Under  this  condition,  the 
test  time  is  controlled  by  the  time  taken  for  the  driver-driven  interface  to 
reach  the  throat,  or  the  leading  expansion  wave  to  deplete  the  reservoir  of 
pressure  behind  the  reflected  shock;  the  flow  duration  is  said  to  be  either 
driver  gas  limited  or  expansion  limited,  respectively.  Figure  80b  shows  the 
flow  duration  in  the  test  section  as  a  function  of  the  Mach  number  of  the 
incident  shock.  Here  it  can  be  seen  that  for  operation  at  low  M^,  running 
times  of  over  25  milliseconds  can  be  obtained  with  a  long  driver  section. 

When  run  under  these  latter  conditions  at  high  pressures  and  Reynolds  numbers. 
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Figure  80b  TEST  TIME  AVAILABLE  FOR  TAILORED-INTERFACE  OPERATION 
OF  48-INCH  SHOCK  TUNNEL 
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the  test  running  times  are  of  the  same  magnitude  or  longer  than  for  piston 
driven  tunnels  (Ref.  15,  16)  with  comparable  stagnation  temperature,  and  the 
reservoir  conditions  and  flow  quality  are  superior  to  piston  driven  tunnels. 
This  is  due  to  the  fact  that  the  test  gas  has  been  processed  by  a  simple  re¬ 
flected  shock  rather  than  multipl  shocks  as  in  piston  tunnels.  A  further  con¬ 
sequence  is  that  the  free  stream  conditions  can  be  calculated  with  far  more 
accuracy  in  a  shock  tunnel . 

It  is  interesting  to  note  that  if  sensitive  high  frequency  instru¬ 
mentation  is  to  be  used  in  the  very  severe  heating  conditions  encountered  in 
turbulent  interaction  regions  in  hypersonic  flow,  running  times  longer  than 
20  milliseconds  present  distinct  problems  because  the  sensing  element  can  be 
damaged  or  destroyed  by  overheating  since  it  must  be  placed  close  to  the  flow 
environment . 

By  running  the  shock  tunnels  at  low  incident  shock  Mach  numbers  and 
high  driver  pressures,  we  can  generate  test  conditions  in  which  we  <.an  obtain 
very  large  Reynolds  numbers.  The  large  test  core  allows  us  to  use  models 
which  are  at  least  3  feet  in  length.  The  maximum  Reynolds  number  capabilities 
of  the  tunnels  are  shown  in  Figure  81 .  Under  these  maximum  Reynolds  number 
conditions,  the  location  of  the  end  of  natural  transition  on  flat  plate 
models  as  determined  from  heat  transfer  and  other  measurements,  is  shown  in 
Figure  82  . 


A. 2  MODEL  AND  FLOW  FIELD  INSTRUMENTATION 

A. 2.1  Introduction 

To  make  a  meaningful  study  of  regions  of  unsteady  flow  over  nose-tip 
configurations,  measurements  of  both  the  mean  and  the  fluctuating  flow  field 
should  be  obtained.  Our  studies  indicated  that  to  obtain  the  complete 
power  spectrum  for  the  surface  pressure,  a  frequency  response  from  200  Hz 
to  100  kHz  was  required.  The  measurements  of  skin  friction  and  heat  transfer 
in  the  separation  and  reattachment  regions  indicated  that  a  frequency  of  at 
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Figure  81  CALSPAN  HYPERSONIC  SHOCK  TUNNEL  PERFORMANCE 
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least  10  kHz  was  required  to  follow  the  motion  of  the  separation  and  reattach¬ 
ment  points.  It  is  also  desirable  to  make  fluctuation  measurements  through 
the  flow  field  to  help  establish  the  source  and  mechanism  whici.  generated  the 
high  level  of  heat  transfer  pressure  fluctuation  observed. 

A. 2. 2  Skin  Friction  Measurements 

In  order  to  define  separation  and  the  length  of  separated  regions 
accurately  in  separated  regions  over  ablated  nose  tip  shapes,  it  is  important 
to  measure  the  surface  shear.  The  inherent  unsteadiness  of  these  regions  also 
makes  it  desirable  to  follow  the  fluctuations  in  the  wall  shear  to  trace  the 
movement  of  the  separation  and  reattachment  points  and  hopefully  gain  some 
insight  into  the  fluid  mechanics  of  such  motions.  Our  recent  studies  indicate 

that  a  frequency  response  of  at  least  10  kHz  is  required  to  follow  these 
movements  in  hypersonic  flow. 

In  axisymraetric  regions  of  viscous  interaction,  separation  is  defined 
as  the  condition  at  which  the  surface  is  zero  at  one  point  only  in  the  inter¬ 
action  region.  Because  direct  measurements  of  surface  shear  are  difficult  to 
make,  pressure  distribution,  surface  pitot,  and  oil  flow  measurements  have 
been  used  to  determine  separation.  Pressure  distribution  measurements  have 
been  shown  by  Green  (Ref.  17)  and  by  Spaid  and  Frishett  (Ref.  18)  to  be  an 
extremely  poor  indicator  of  boundary  layer  separation  in  supersonic  flow.  Holden 
(Ref.  19)  found  this  technique  was  even  more  insensitive  at  hypersonic  speeds. 
From  our  studies  we  have  noted  that  t'  e  introduction  of  a  small  disturbance  into 
the  laminar  sublayer  beneath  the  sonic  line  can  cause  dramatic  effects  both 
in  the  immediate  vicinity  and  downstream  of  the  disturbance. 

On  our  studies  of  the  separation  of  turbulent  boundary  layers  in 
hypersonic  flow,  we  have  used  a  gage  (Ref.  20)  which  directly  measures  surface 
shear  to  indicate  the  separation  condition.  A  diagram  of  the  skin  friction 
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transducer  which  we  would  use  in  the  proposed  studies  is  shown  in  Figure  83. 
The  transducer  consists  of  a  diaphragm  which  is  supported  flush  with  the  model 
surface  by  two  piezo-ceramic  beams,  which  develop  a  charge  when  placed  in 
bending  by  a  surface  shear  on  the  diaphragm.  A  third  beam  is  used  to  provide 
acceleration  compensation;  the  beams  are  connected  electrically  to  eliminate 
thermal,  normal,  and  transverse  pressure  effects.  An  FET  impedance  transform 
circuit  IS  mounted  internally  to  eliminate  cable  noise  effects  at  low  levels 
of  skin  friction.  The  gage,  which  has  been  refined  and  developed  over  the 
past  12  years,  has  been  used  to  measure  very  low  levels  of  skin  friction 
encountered  in  separated  regions  in  low  Reynolds  number  hypersonic  flow  and 
more  recently  very  high  levels  in  regions  of  shock  wave-turbulent  boundary 
layer  interactions  in  hypersonic  flow.  Because  of  the  very  severe  heating 
conditions  encountered  in  the  latter  studies,  special  care  was  taken  to 
minimize  the  heat  conduction  through  the  flexures.  The  very  large  dynamic 
loads  generated  on  the  transducers  during  tunnel  shutdown  when  run  at  the 

high  dynamic  pressure  conditions  used  in  our  studies  caused  the  diaphragms  to 
be  torn  from  the  supporting  beam.  This  problem  was  overcome  by  careful  design 
of  the  flexure  and  by  mounting  the  transducer  in  the  seismic  mass-rubber 
suspension  system  shown  in  Figure  83. 


A. 2. 3  Heat  Transfer  Instrumentation 


A  knowledge  of  the  heat  transfer  distribution  on  nose  tips  is  of 
great  importance  because  of  the  very  severe  heating  rates  generated  in  the 
reattachment  region.  Almost  as  important  as  the  severity  of  reattachiiient 
heating  is  the  extremely  large  heat  transfer  gradients  which  occur  both  in 
the  separation  and  reattachment  regions.  Regions  of  high  heat  transfer 
gradient  present  a  problem  to  the  experimenter  because  they  can  cause  trans¬ 
verse  heat  conduction  problems  in  the  model.  This  can  distort  heat  transfer 
distributions  and  in  some  cases  dramatically  reduce  the  maximum  indicated 
heating.  An  example  which  illustrates  this  feature  is  given  in  the  report 
of  Hiers  and  Loubsky  (Ref.  21)  where  they  used  a  thin  skin  technique  to  obtain 
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heat  transfer  measurements  in  regions  of  turbulent  viscous  layer  attachment. 

They  found  that  the  corrections  to  their  measurements  were  of  the  same 
magnitude  as  the  measurements  themselves.  To  reduce  transverse  heat  conduc¬ 
tion  along  the  model  surface,  the  model  should  be  made  from  a  non-conducting 
material;  this  is  incompatible  with  the  thin  skin  technique  unless  a  segmented 
model  is  used.  To  overcome  these  problems  we  use  a  measurement  of  heat 
transfer  which  relies  on  sensing  the  transient  surface  temperature  of  a  non¬ 
conducting  model  by  means  of  thin-film  resistance  thermometers.  Because  the 
thermal  capacity  of  the  gage  is  negligible,  the  instantaneous  surface  tempera¬ 
ture  of  the  backing  material  is  related  to  the  heat  transfer  rate  by  the 
classical  semi-infinite  slab  theory.  Analog  networks  are  used  to  convert  the 
outputs  of  the  gages,  which  are  proportional  to  surface  temperature,  to  a 
voltage  directly  proportional  to  heat  transfer.  The  thin  film  gage  has  a 
frequency  response  to  1  MHz. 

The  gages  are  fabricated  on  either  small  pyrex  buttons  or  on  contoured 
inserts,  the  ability  of  this  technique  to  make  closely  spaced  measurements  in 
regions  of  shock  impingement  is  demonstrated  in  Figure  81. 

A. 2. 4  Surface  and  Flow  Field  Pressure  Measurements 


We  employed  two  types  of  surface  pressure  transducers  in  our  shock 
tunnel  studies.  The  Calspan-designed  and  constructed  lead  zirconium  titanate 
piezoelectric  pressure  transducers  (Ref-  22)  were  used  to  obtain  essentially 
the  mean  pressure  distribution  through  the  interaction  region,  though  the 
transducer  and  orifice  combination  could  follow  fluctuations  up  to  15  kHz. 

A  second  flush-mounted  transducer,  especially  designed  for  high  frequency 
measurements  by  PCB  in  Buffalo,  was  used  to  obtain  surface  pressure  fluctua¬ 
tion  measurements  from  200  Hz  to  120  kHz.  To  prevent  a  resonance,  a  special 
mounting  system  was  developed  (as  shown  in  Figure  84)  to  lock  the  gage  firmly 
into  the  model.  A  thin  insulating  barrier  of  aluminized  mylar  was  attached 

to  the  diaphragm  of  the  transducer  to  prevent  thermal  heating  effect  caused 

2 

by  the  large  heat  transfer  rates  (over  500  Btu/ft  sec)  generated  in  the 
reattachment  regions . 
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Figure  83  DRAWING  OF  SECTION  THROUGH  SKIN  FRICTION  TRANSDUCER 


heat  shield  (ALUMIHIZED  mylar  .002  THICK) 


(a)  TYPICAL  HOUHTIHG  TECHHIQUE  USED  IN  MODEL 
Figure  84  HIGH-FREQUENCY  PRESSURE  MOUNTING 
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Figure  85  MODELS  SHOWING  TYPICAL  HEAT  TRANSFER  GAGE  INSTALLATION 
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A.  2.5 


Data  Recording  and  Processing; 


The  outputs  from  tJie  transducers  are  recorded  on  a  NAVCOR  (Ref.  14j 
magnetic  drum  system  and  on  CliC  and  AMRRX  DM  tape  recorders,  and  monitored 
on  oscilloscopes.  The  NAVfOR  system,  which  can  hold  48  channels  of  data  in 
digital  form,  is  essentially  a  low  frequency  response  system  which  will  record 
the  mean  output  of  the  gages  over  the  running  time.  The  fluctuation  measure¬ 
ments  will  he  recorded  on  the  FM  tape  recorders  in  analog  form  and  subsequently 
converted  into  digital  form  by  the  analog-to-digital  conversion/data  storage 
system  developed  at  falspan. 

Ihe  fluctuation  data,  in  punched  card  form,  are  then  processed  by 
an  existing  digital  computer  program  to  obtain  the  statistical  properties  -- 
Root  Mean  Square  Pressure,  Power-  and  Cross-Spectra,  and  Auto-  and  Cross- 
Correlation  functions. 

A. 2.6  Single  Frame  and  High  Speed  Motion  Schlieren  Photography 

A  single-pass  Schlieren  system  with  a  horizontal  knife-edge  was  em¬ 
ployed  in  all  of  the  studies  of  high  Reynolds  number,  high  Mach  number  flows 
over  ablated  nose  shapes.  However  for  the  studies  of  the  low  Reynolds  numbers 
aerodynamics  of  the  Widhopf  1  and  2  nose  tips  we  employed  a  double-pass  sys¬ 
tem  with  a  high  degree  of  cut-off.  Here  we  used  a  single  frame  camera  with 
a  microsecond  spark  to  record  the  event.  We  used  a  Red  Lake  high  speed  movie 
camera  with  speeds  of  up  to  10,000  frames  per  second  to  record  the  flow  oscil¬ 
lation  occurring  about  the  highly  ablated  nose  shapes  and  the  particle- launch 
model. 
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CONCLUSIONS 


A  unique  set  of  experimental  measurements  have  been  reported  which 
demonstrate  that  a  single  minute  particle  launched  through  the  bow  shock  of  a 
flat-ended  cylinder  can  augment  the  heating  rates  to  the  model  surface  by  as 
much  as  a  factor  of  ten.  These  studies  have  identified  the  principal  mech¬ 
anism  of  enhancement  heating  as  one  associated  with  the  reattachment  of  a 
jet/shear  layer  generated  by  the  shock-shock  interaction  after  the  particle 
has  penetrated  the  bow  shock.  This  investigation  identified  four  distinctly 
different  classes  of  flow  phenomena  resulting  from  particle  shock-bow  shock 
interaction,  together  with  the  heating  augmentation  levels  associated  with 
them.  The  most  dramatic  of  these  is  one  where  a  flow  instability  similar  to 
that  encountered  over  the  highly  indented  nose  shapes  is  observed. 
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